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ABSTRACT 

An  analytical  study  of  icing  simulation  for  turbine  engines  in  altitude  test  cells  was 
made.  A  mathematical  model  of  a  typical  direct-connect  type  of  icing  test  cell  was 
developed,  and  the  governing  equations  were  programmed  for  computer  solution  in 
FORTRAN  computer  language.  A  parametric  study  was  performed  to  determine  the  effects 
of  the  test  cell  inlet  and  water  spray  conditions  on  the  thermodynamic  and  kinetic  state 
of  the  flow  in  the  test  section  or  at  the  engine  compressor  face.  The  importance  of  correctly 
simulating  droplet  size  distribution,  mean  effective  droplet  diameter,  liquid  water  content, 
and  humidity  was  investigated.  The  results  of  the  investigation  lend  further  support  to 
the  fact  that  ground  test  facilities  provide  the  best  available  capability  for  conducting 
turbine  engine  icing  tests.  The  ability  to  define  and  control  the  simulated  environment 
gives  an  altitude  test  cell  distinct  advantages  not  enjoyed  by  flight  testing  in  natural  icing 
environments  or  flight  testing  in  artificial  environments  created  by  tanker  aircraft. 
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SECTION  I 
INTRODUCTION 


1.1  THE  ICING  PROBLEM 

Ice  accretion  on  aircraft  usually  causes  a  degradation  of  the  aircraft  performance 
and  operational  safety.  The  degree  of  degradation  is  related  to  the  amount  of  ice  adhering 
to  the  aircraft  and  propulsion  system  components.  Because  of  this,  aircraft  are  routinely 
equipped  with  either  a  system  for  intermittently  removing  the  accumulated  ice  (a  de-icing 
system)  or  a  system  for  continuously  maintaining  the  exposed  surfaces  of  the  aircraft 
free  from  ice  (an  anti-icing  system).  These  systems  must  be  tested  under  conditions  which 
closely  simulate  or  if  possible  duplicate  the  natural  atmospheric  icing  conditions.  This 
requires  that  natural  icing  conditions  be  well  known  and  that  methods  be  developed  for 
simulating  these  conditions  for  testing  purposes.  It  follows  that,  to  ensure  good  simulation, 
the  primary  variables  which  define  icing  conditions  must  be  established  and  controlled 
in  the  test  facility  during  the  testing  process. 

A  propulsion  engine  test  facility  (Fig.  1,  Appendix  I)  provides  a  unique  capability 
for  simulating  aircraft  engine  environments  since  the  aerodynamic  and  thermodynamic 
conditions  at  the  compressor  face  are  accurately  controllable.  In  this  facility,  the  engine 
is  coupled  directly  to  the  air  supply  ducting,  and  the  attempt  is  made  to  simulate  the 
flow  conditions  that  would  exist  at  the  compressor  face  in  flight.  This  type  of  facility 
has  been  utilized  for  conducting  icing  tests ;  of  turbine  engines.  A  recent  survey  of  the 
state-of-the-art  of  icing  simulation  for  propulsion  engine  testing  was  reported  by  Kissling 
(Ref.  1).  An  important  conclusion  of  Kissling’s  is  that  the  ground  test  facility  offers  the 
best  solution  for  determing  anti-icing  system  performance  because  all  variables  which  define 
an  icing  condition  can  be  very  accurately  controlled. 

The  purpose  of  this  study  was  to  analyze  the  simulation  of  icing  conditions  in 
propulsion  engine  test  cells  as  currently  practiced.  Only  the  direct-connect  type  of  test 
cell  was  considered.  While  this  limits  the  applicability  of  the  results  of  the  analysis  primarily 
to  the  testing  of  pod-mounted  engines,  the  basic  approach  is  applicable  to  fuselage-mounted 
engines  as  well.  In  the  case  of  engines  having  inlet  ducts,  as  is  usually  the  case  with 
fuselage-mounted  engines,  it  is,  of  course,  necessary  to  consider  the  effect  of  the  inlet 
ducting  on  the  flow  conditions  at  the  compressor  face. 

1J2  TOPICS  CONSIDERED  IN  THE  REPORT 

Section  II  of  this  report  includes  discussions  of  natural  icing  conditions  and  the 
military  specifications  for  conducting  icing  tests  for  turbine  engines.  In  Section  III,  some 
of  the  problems  associated  with  simulating  icing  conditions  in  an  altitude  test  cell  are 
discussed.  An  analysis  of  flow  in  a  direct-connect  type  of  test  cell  is  presented  in  Section 
IV,  and  an  application  of  the  analysis  follows  in  Section  V.  In  Section  VI,  a  number 
of  factors  affecting  icing  simulation  are  discussed. 
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SECTION  II 

DISCUSSION  OF  NATURAL  ICING  CONDITIONS 
2.1  METEOROLOGICAL  DATA 

The  principal  factors  which  determine  an  icing  condition  are  cloud  temperature,  liquid 
water  content,  and  water  droplet  size  distribution,  including  the  "mean  effective"  droplet 
diameter  for  the  given  droplet  size  distribution.  In  addition,  the  extent  and  severity  of 
the  cloud  formation  causing  the  icing  condition  are  important  factors,  i.e.,  cloud  size  and 
shape.  The  NACA  investigated  atmospheric  icing  phenomena  and  reported  the  data  from 
these  studies  during  the  period  of  years  between  1940  and  1960.  References  2  through 
5  are  typical  reports  from  this  era  of  icing  investigation.  Flights  were  made  with  specially 
instrumented  aircraft  through  weather  conditions  that  led  to  aircraft  icing.  The  data  from 
these  flights  were  cloud  liquid  water  content,  cloud  temperature  as  a  function  of  the 
pressure  altitude,  and  the  relation  between  cloud  formation  or  configuration  and  modes 
of  ice  accretion.  The  modes  of  ice  accumulation  were  defined  by  the  rates  of  ice  buildup 
on  the  aircraft,  that  is,  if  the  rates  were  small  or  large,  and  if  the  accumulation  was 
continuous  or  intermittent  in  time.  Also,  the  shape  and  the  texture  of  the  ice  accumulation 
were  studied. 

Essentially,  it  was  found  that  the  severity  of  icing  can  be  correlated  with  the  type 
of  icing  cloud  formation.  Moreover,  intermittent  icing  can  be  associated  with  one  type 
of  cloud  formation,  and  continuous  icing  can  be  associated  with  another  type  cloud 
formation.  For  clarity,  the  method  of  summarizing  data  which  was  presented  in  Ref.  4 
is  quoted  in  this  report  as  follows:  In  Ref.  4,  it  is  stated  that  "the  data  on  liquid-water 
content,  mean-effective  droplet  diameter,  ambient  air  temperature,  and  pressure  altitude 
of  icing  conditions  are  divided  into  two  categories:  for  stratiform  and  cumuliform  clouds. 
Stratiform  clouds  are  layer-type  clouds  in  which  the  horizontal  extent  is  greater  than  the 
vertical  thickness.  This  cloud  type  includes  stratus,  strato-cumulus,  alto-stratus,  and 
nimbo-stratus.  Cumuliform  clouds  are  clouds  with  vertical  extent  comparable  to  the 
horizontal  extent.  This  cloud  type  includes  cumulus,  cumulus  congestus,  cumulo-nimbus, 
and  alto-cumulus-castellatus.  This  division  of  icing  situations  was  necessary  because  the 
icing  conditions  in  the  two  cloud  types  are  very  different,  owing  to  the  difference  in 
the  manner  in  which  they  are  formed."  The  continuous  mode  of  icing  is  associated  with 
the  stratiform  clouds,  and  the  intermittent  icing  mode  is  associated  with  the  cumuliform 
cloud  type.  On  the  average,  the  liquid  water  content  (grams  of  water  per  cubic  meter 
of  cloud)  was  found  to  be  greater  in  the  cumuliform  cloud  than  in  the  stratiform  cloud 
formations.  Thus,  moderate  to  heavy  but  short  duration  icing  was  encountered  on  the 
average  in  cumuliform  clouds  and  light  to  moderate  but  long  duration  icing  was  incurred 
on  the  average  in  the  stratiform  cloud  types.  A  concise  summary  of  the  meteorological 
data  for  icing  conditions  which  occurs  on  the  average  for  the  cumuliform  and  the  stratiform 
clouds  is  provided  in  Ref.  6  and  is  reproduced  and  presented  in  Figs.  2  through  5.  Figures 
2  and  4  present  the  correlation  between  mean  effective  droplet  diameter  and  liquid  water 
content  for  stratiform  and  cumuliform  clouds,  respectively.  Figures  3  and  5  present  the 
relation  between  ambient  temperature  and  pressure  for  stratiform  and  cumuliform  clouds, 
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respectively.  It  should  be  noted  that  the  iqing  conditions  defined  at  the  time  the  data 
were  being  obtained  were  defined  in  terms  of  typical  aircraft  speed,  rate  of  climb,  and 
performance  in  that  era.  Icing  conditions  for  modem  aircraft  may  be  redefined  in  the 
light  of  modern  aircraft  performance.  That  is,  modem  aircraft  performance  may  well  be 
sufficient  reason  to  redefine  continuous  or  intermittent  icing  and  the  relative  severity 
associated  with  each  type  icing  condition.  It  is  also  important  to  note  that  the  FAA 
specifications  (Ref.  6)  for  withstanding  icing  pertain  specifically  to  icing  on  external 
aerodynamic  surfaces  and  not  to  engine  performance  under  icing  conditions. 

Despite  the  fact  that  there  is  still  much  of  cloud  physics  which  remains  to  be 
investigated,  sufficient  data  have  been  obtained  for  the  icing  problem  to  enable  the 
construction  of  a  fairly  comprehensive  description  of  the  icing  environment.  According 
to  Lewis  (Ref.  7),  the  amount  of  data  gathered  from  reported  icing  conditions  covers 
such  a  range  and  is  sufficiently  complete  that  further  investigation  of  natural  environments 
is  unnecessary.  If  this  is  accepted  then  the  data  in  the  FAA  circular  can  be  considered 
as  adequately  describing  the  envelope  of  possible  icing  environments  for  all  design  and 
test  purposes. 

2.5  MILITARY  SPECIFICATIONS  FOR  TURBOJET  AND  TURBOFAN  ENGINE 
TESTING  UNDER  ICING  CONDITIONS 

From  the  NACA  data,  such  as  summarized  by  the  FAA  circular,  standard  icing 
conditions  have  been  defined  for  the  purpose  of  testing  turbojet  and  turbofan  engines 
for  military  applications.  These  conditions  are  specified  in  Ref.  8  and  are  presented  in 
Table  I  (Appendix  II).  The  military  specifications  appear  to  be  definitive  at  first  glance. 
They  define  the  required  gas  temperature,  the  liquid  water  content,  the  flight  Mach  number, 
and  the  geometric  altitude  from  which  the  pressure  altitude  can  be  obtained.  Also  specified 
is  the  "mean  effective"  diameter  of  the  liquid  droplets  in  the  icing  cloud.  Unfortunately, 
there  is  ambiguity  associated  with  these  specifications.  For  example,  the  definition  of  the 
"mean  effective"  drop  diameter  is  not  precisely  clear.  There  have  been  two  commonly 
defined  mean  effective  drop  diameters.  One  is  defined  as  the  mass  or  volume  mean  diameter. 
The  volume  mean  diameter  is  defined  by  the  equation: 

/Nl 

d  =  (E  (Nj/TVdj3 

where  Nj  equals  the  number  of  drops  having  a  particular  diameter  (dj)  in  a  given  sample 
and  Nt  is  the  sum  total  number  of  all  drops  in  the  sample.  On  the  other  hand,  another 
mean  effective  diameter  is  the  mass  median  diameter  defined  as  the  diameter  (dm )  of 
those  droplets  for  which  one-half  of  the  entire  liquid  water  content  in  the  spray  cloud 
is  contained  in  droplets  whose  diameter  is  less  than  dm.  The  values  taken  by  d  and  dm 
from  a  given  sample  depend  on  the  distribution  of  the  mass  of  liquid  water  among  the 
droplets  in  accordance  to  their  diameters.  Thus,  the  statement  that  the  mean  effective 
diameter  should  be  "M"  microns  for  icing  testing  is  ambiguous  unless  the  definition  of 
the  mean  effective  diameter  is  also  given.  It  should  be  noted  that  most  NACA  publications 
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and  the  major  engine  manufacturers  use  the  mass  median  diameter  (dm)  as  the  mean 
effective  diameter.  It  may  be  necessary  to  specify  the  droplet-size  distribution  for  icing 
tests  where  the  amount  of  ice  accumulated  is  a  strong  function  of  the  droplet  size.  This 
effect  will  be  discussed  in  a  following  section  of  the  report.  In  this  report,  the  mass  median 
diameter  is  used  as  the  mean  effective  diameter. 

Another  possible  source  of  ambiguity  is  in  the  definition  of  gas  temperature.  The 
gas  temperature  is  specified  as  both  an  atmospheric  air  temperature  and  as  an  "inlet" 
temperature.  The  terminology  is  unclear,  since  engines  are  tested  with  and  without  an 
inlet  duct.  In  flight,  an  engine  which  is  pod-mounted  is  exposed  directly  to  the  free  stream. 
Hence,  the  direct -connect  testing  of  such  an  engine  is  conducted  by  ducting  air  to  the 
compressor  face,  which  is  conditioned  to  simulate  conditions  experienced  at  the  compressor 
face  during  flight.  A  turbojet  engine  for  a  fighter-type  aircraft  typically  has  the  engine 
mounted  in  the  interior  of  the  aircraft  (fuselage-mounted),  and  the  free-stream  airflow 
is  ducted  to  the  engine  from  some  inlet  air  scoop.  Thus,  there  are  inlet  conditions  to 
the  inlet  ducting  and  inlet  conditions  to  the  engine  compressor.  Therefore,  the  term  inlet 
temperature  in  the  military  specifications  is  not  sufficiently  clear.  A  burden  is  put  on 
the  test  engineer  to  decide  which  inlet  condition  to  consider.  In  previous  icing  tests,  it 
has  been  assumed  that  the  specified  temperature  is  the  free-stream  total  temperature.  The 
ambiguity  in  the  military  specification  can  be  removed  if  the  defined  conditions  are 
specified  as  free-stream  total  conditions  rather  than  inlet  conditions.  Thus  T„  or  Tt  p„ 
or  p(^,  humidity  (co„,),  liquid  water  content  (L),  M„„,  and  d  or  dm  can  completely  define 
icing  conditions.  Note  that  in  this  study  it  is  recommended  that  the  free-stream  humidity 
(cj„)  be  also  specified.  Reasons  for  this  recommendation  will  be  given  in  the  following 
sections  of  this  report.  At  present,  u>„  is  not  given  in  the  military  specification.  In  typical 
previous  tests,  has  been  set  to  the  saturation  value  corresponding  to  Tt^  and  pt>o.  It 
is  expected  that,  in  natural  icing  conditions,  the  relative  humidity  is  approximately  100 
percent  so  that  cu„  corresponds  to  saturation  at  the  ambient  atmospheric  conditions. 

Another  important  point  to  consider  for  defining  atmospheric  icing  conditions  is 
whether  the  drop-size  distribution  must  be  specified  and  required  during  test  simulation 
in  addition  to  specification  of  the  mean-effective  diameter.  The  reason  is  that  the  icing 
on  engine  components  could  be  sensitive  to  drop-size  distribution.  This  question  is 
considered  later  in  the  report. 

To  summarize  this  section  of  the  report,  the  military  specifications  for  icing  testing 
can  be  made  sufficiently  definitive  by  specifying  free-stream  static  or  total  conditions  of 
pressure  and  temperature,  flight  Mach  number,  free-stream  humidity,  liquid  water  content, 
and  a  precisely  defined  mean  effective  diameter.  This  latter  parameter  should  be  specified 
as  either  the  mean  volumetric  or  mass  median. 
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SECTION  III 

ICING  CONDITIONS  IN  PROPULSION  ENGINE  TEST  CELLS 
3.1  DISCUSSION 

Ideally,  a  test  cell  for  conducting  icing  tests  should  duplicate  the  flow  conditions 
at  the  compressor  face  experienced  by  an  engine  in  flight  through  an  icing  cloud.  The 
flow  conditions  are  characterized  by  temperature,  humidity,  pressure,  velocity,  liquid  water 
content,  mean  effective  droplet  size,  and  droplet-size  distribution.  It  is  generally  not 
practical  and  frequently  impossible  to  duplicate  all  of  the  flow  conditions  encountered 
in  atmospheric  icing.  Fortunately,  however,  it  is  not  necessary  to  duplicate  all  of  the 
icing  conditions  to  get  adequate  simulation. 

Briefly,  the  method  used  to  produce  ice  accretion  on  turbojet  and  turbofan  engine 
inlet  components  for  icing  testing  is  to  inject  a  continuous  spray  of  water  droplets  into 
a  cold  airstream  directed  at  the  engine.  A  typical  icing  test  cell  is  shown  schematically 
in  Fig.  1.  The  water  droplets  are  accelerated  to  near  the  airspeed  by  aerodynamic  drag. 
Through  heat  and  mass  transfer,  the  droplets  are  expected  to  come  to  the  proper 
supercooled  condition.  Air-atomizing  spray  nozzles  are  typically  used  to  produce  the  spray 
cloud  of  droplets  having  the  desired  mean  effective  diameter. 

In  the  previous  section,  it  was  noted  that  the  icing  specifications  are  given  in  terms 
of  flight  conditions.  Since  the  flow  conditions  at  the  compressor  face  are  rarely  the  same 
as  the  flight  conditions,  the  air  having  been  accelerated  or  diffused  before  entering  the 
engine,  it  is  necessary  to  relate  the  engine  inlet  conditions  to  the  free-stream  conditions. 
This  is  discussed  further  in  the  following  sections.  The  problem  of  establishing  and 
monitoring  the  liquid  water  content  and  mean  effective  droplet  size  of  the  spray  cloud 
in  the  test  cell  is  considered  in  the  following  section. 

3 2  RELATION  BETWEEN  FREE-STREAM  AND  ENGINE 

COMPRESSOR  INLET  CONDITIONS 

For  purposes  of  discussion,  flow  from  the  free  stream  to  the  engine  compressor  inlet 
is  assumed  to  be  one-dimensional  and  isentropic.  This  is  a  good  assumption  for  a 
pod-mounted  engine.  The  validity  of  the  assumption  for  an  engine  having  a  long  inlet 
duct  depends  on  the  duct  shape  and  flow  conditions. 

To  simulate  the  free-stream  temperature,  pressure,  and  flight  velocity  in  the  test  cell, 
the  total  pressure  and  total  temperature  at  the  compressor  inlet  are  set  equal  the  total 
pressure  and  total  temperature  computed  from  the  altitude  static  pressure  and  static 
temperature  and  the  flight  Mach  number.  Based  on  results  to  be  discussed  later,  it  has 
been  found  that  flow  from  the  free  stream  to  the  compressor  inlet  occurs  at  essentially 
constant  specific  humidity.  Hence,  the  specific  humidity  at  the  compressor  inlet  in  the 
test  cell  should  be  set  equal  to  the  specific  humidity  of  the  air  at  the  flight  conditions 
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being  simulated.  In  a  natural  icing  cloud,  it  is  expected  that  the  supercooled  liquid  is 
in  equilibrium  with  its  saturated  vapor  at  the  local  ambient  temperature.  That  is,  the 
relative  humidity  in  the  cloud  is  100  percent.  The  specific  humidity  (go)  at  the  compressor 
is  given  by 


co  =  = - =  0.622  - 

Ma  P«  P« 

where  pVoo  is  the  saturation  pressure  of  supercooled  liquid  at  the  free-stream  static 
temperature  (T„),  and  p„  is  the  free-stream  static  pressure,  i.e.,  altitude  pressure.  The 
relative  humidity  at  the  compressor  will  be  1  00  percent  only  if  the  compressor  inlet  flow 
velocity  equals  the  flight  speed.  Thus,  the  flow  entering  the  engine  will  generally  not 
be  in  thermal  equilibrium  since  equilibrium  can  exist  only  at  100-percent  relative  humidity. 
The  importance  of  humidity  on  the  icing  problem  is  that  humidity  affects  the  heat  and 
mass  transfer  from  droplets  in  the  spray  and  the  heat  and  mass  transfer  from  exposed 
surfaces  within  the  engine  inlet.  These  effects  are  analyzed  and  illustrated  quantitatively 
in  the  following  sections. 

Acceleration  or  deceleration  of  the  flow  from  the  free  stream  to  the  compressor 
inlet  results  in  a  change  in  the  liquid  water  content.  Two  separate  phenomena  are  involved. 
One  is  the  compressibility  of  the  gas.  The  second  phenomenon  is  droplet  inertia  and  the 
inability  of  the  droplets  to  follow  streamlines  exactly.  This  problem  is  discussed  in  detail 
in  a  following  section.  It  is  sufficient  to  remark  here  that  the  liquid  water  content  generally 
changes  by  less  than  30  percent  from  free  stream  to  engine  face.  In  Section  VI,  it  is 
shown  that  the  power  required  to  maintain  an  anti-icing  condition  is  not  highly  dependent 
on  liquid  water  content.  Thus  a  30-percent  correction  to  the  liquid  water  content  may 
not  be  an  important  consideration  for  icing  cloud  simulation. 

3.3  MEAN  EFFECTIVE  DROPLET  DIAMETER  DETERMINATION 

During  engine  icing  tests  in  propulsion  facilities,  it  is  necessary  to  ensure  that  the 
water  spray  cloud  is  characterized  by  the  desired  mean  effective  drop  diameter.  One 
established  method  to  measure  drop  sizes  is  to  capture  a  sample  of  drops  from  the  spray 
cloud  on  an  oil-film  covered  slide.  The  slide  is  usually  exposed  to  the  water-air  mixture 
for  a  short  time  at  some  flow  station  downstream  of  the  water  spray  station.  Then  the 
slide  water  droplet  sample  is  analyzed  to  determine  the  droplet  size-frequency  distribution. 
From  this  distribution,  the  mean  effective  diameter  can  be  calculated.  Unfortunately,  this 
method,  while  simple  and  convenient,  is  not  accurate.  Correction  factors  must  be 
introduced  into  the  analysis  of  the  droplet  distribution  because  of  droplet  evaporation 
after  being  captured,  because  the  drops  may  not  be  spherical  in  form  in  the  oil  film 
but  may  flatten,  and  because  the  slide  samples  only  a  very  small  portion  of  the  air-water 
mixture  flowing  through  the  duct.  A  much  more  accurate  and  modem  technique  for 
determining  conditions  in  the  water  spray  cloud  utilizes  laser  holography.  With  a 
laser-holographic  camera  system  such  as  described  by  Gall  and  Floyd  (Ref.  9)  a 
three-dimensional  photograph  of  any  point  in  the  duct  flow  can  be  obtained.  From  the 
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photograph,  an  analysis  of  the  number  of  drops,  their  size,  and  free-flight  shape  can  be 
made.  With  a  sufficiently  varied  sampling  across  the  entire  duct,  an  accurate  analysis  of 
the  spray  cloud  can  be  obtained.  This  process  can  be  done  "on-line"  during  the  icing 
test,  using  an  electronic  scanning  device  to  count  only  those  water  droplets  whose  form 
is  clearly  in  focus  in  the  hologram. 

a4  TECHNIQUE  FOR  SETTING  LIQUID  WATER  CONTENT 


The  proper  liquid  water  content  of  the  air-water  mixture  is  maintained  by  metering 
the  airflow  and  the  water  flow,  and  adjusting  the  water  flow  rate  to  the  desired  amount. 
The  relationship  between  the  liquid  water  content  (L  or  £)  and  the  actual  water  loading 
factor  (fc)  in  lbm  of  water  per  Ibm  of  air  can  be  obtained  from  considering  the  flow 
in  a  stream  tube  in  the  free  stream. 
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or  in  particular  where  the  liquid  water  content  is  in  gm/m3,  T  in  °R,  and  p  in  psf: 

L  T 

f  =  3.326  x  10  3  — - -  (2) 

°oo  P.o 

Since  L  is  usually  specified  for  free-stream  icing  conditions,  T  and  p  are  the  free-stream 
static  temperature  and  pressure,  respectively.  The  engine  mass  flow  in  the  test  facility 
is  metered,  often  using  a  critical-flow  venturi.  In  this  case, 


W 


a 


0.532  Cj  f  p(  A ( 

*  ©O 


oo 


(3) 


for  mass  flow  metering  venturis  operating  at  a  critical  (choked)  condition;  Qj  represents 
a  known  discharge  or  flow  coefficient.  By  using  this  equation  and  Eq.  (2)  for  fc,  the 
water  flow  rate  can  be  written  as 

0.532  Cd  fpt  A ( 

W  =  -  -°°- - 
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Since  T„  and  p„  are  related  to  the  total  conditions  (Tt^  and  ptoJ  through  the  isentropic 
Mach  number  relations,  the  equation  for  Ww  can  be  written: 
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where  M„  is  the  free-stream  flight  Mach  number  for  the  icing  condition  being  simulated. 
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Equation  (1)  is  valid  so  long  as  the  water  in  the  streamtube  remains  constant.  In 
flight,  the  engine  often  causes  the  flow  upstream  of  the  inlet  to  accelerate  or  decelerate 
depending  on  the  particular  flight  condition  (see  Fig.  6).  If  the  water  droplets  have 
sufficient  inertia,  it  is  possible  that  they  would  not  follow  streamlines  entering  the  engine. 
For  this  case,  Eq.  (1)  does  not  give  the  correct  loading  factor.  The  amount  by  which 
Eq.  (1)  would  have  to  be  corrected  when  the  droplet  inertia  is  important  can  be 
quantitatively  estimated  as  follows:  The  basic  equation  for  the  loading  factor  was 


The  amount  of  water  (Ww)  to  be  sprayed  in  the  test  cell  is  related  to  the  conditions 
in  the  atmosphere  encountered  by  the  engine  during  its  operation.  Clearly,  Ww  thus  depends 
on  the  free-stream  liquid  water  content  (£„)>  the  flight  velocity,  and  whether  droplets 
follow  streamlines  or  not.  When  the  engine  is  operating  under  static  conditions,  it  is  quite 
reasonable  to  expect  that  the  amount  of  water  ingested  is  equal  to  the  liquid  water  content 
(£_)  times  the  air  volumetric  flow  rate  through  the  engine.  Hence, 
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Therefore,  it  is  not  necessary  to  consider  droplet  inertia  for  the  static  ambient  icing 
conditions,  and  this  case  is  identical  to  the  case  where  droplets  follow  streamlines. 


The  other  case  to  consider  is  where  droplets  follow  straight  paths,  i.e.,  their  trajectories 
are  not  influenced  by  the  disturbed  flow  field  upstream  of  the  engine  inlet.  In  this  case, 
the  water  flow  rate  is  given  by 
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On  the  other  hand,  the  airflow  rate  is  given  by 
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But  since  fc  can  be  defined  by 
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as  the  free-stream  loading  factor,  then 
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(12) 


This  is  the  required  loading  factor  when  droplets  are  not  influenced  by  streamlines.  To 
put  this  in  a  more  precise  form,  the  flow  from  the  free  stream  to  the  compressor  face 
is  approximately  isentropic,  in  which  case 


Thus, 
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This  relation  has  been  evaluated  for  the  military  icing  specifications,  assuming  two 
representative  values  of  Mcf,  and  the  results  are  presented  in  Fig.  7.  The  relationship  is 
shown  between  the  water  loading  factor  which  must  be  established  in  the  test  cell  to 
the  free-stream  water  loading  factor  based  on  free-stream  liquid  water  content,  pressure, 
and  temperature.  The  two  curves  represent  limiting  cases  for  establishing  the  water  loading 
factor  for  icing  tests.  Where  droplets  follow  streamlines  exactly  and  for  static  environment 
icing  conditions,  fc  =  fCoc.  Where  droplets  are  not  influenced  by  streamlines,  fc  =  T  fc_, 
where  T  must  be  determined  by  the  values  of  and  Mcf. 


A  question  arises  as  to  which  limiting  case  is  more  realistic  for  evaluating  fc  for 
testing  purposes;  that  is,  do  droplets  follow  streamlines  or  not?  To  estimate  representative 
values  of  T,  an  average  T  value  can  be  defined  as 


r  = 


r+  i 

2 


(15) 


and  the  following  table  was  devised: 


MCf 

r  =  fc/fc. 

r 

0.32 

0.25 

1.25 

1.125 

0.71 

0.25 

2.21 

1.605 

0.32 

0.50 

0.70 

0.85 

0.71 

0.50 

1.23 

1.1 15 
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It  is  evident  from  this  table  that,  on  the  average  for  military  testing  purposes,  the  highest 
amount  of  additional  water  which  must  be  injected  into  a  test  cell  to  simulate  correctly 
the  proper  engine  ingestion  of  water  is  about  60  percent  greater  than  the  free-stream 
water  loading.  This  pertains  to  the  case  where  the  engine  is  throttled  back,  yet  the  airspeed 
remains  fairly  high  during  aircraft  descent.  This  result  compares  with  the  "scoop  factor" 
of  1.33  presented  in  Ref.  10  for  the  same  kind  of  flight  condition  called  the  "idle  descent" 
mode.  For  the  majority  of  the  test  points,  however,  as  indicated  by  the  average  T  factor, 
the  test  cell  loading  factor  (fc)  is  approximately  equal  to  the  free-stream  loading  factor. 
The  next  important  question  is  to  determine  how  sensitive  the  engine  anti-icing  power 
requirements  are  to  variations  in  liquid  water  content,  hence  to  fc.  If  they  are  relatively 
insensitive,  then  it  is  reasonable  to  establish  fc  =  fc_  (Eq.  (2))  for  all  icing  tests  defined 
in  the  military  specifications. 

3.5  SUMMARY 

The  air  at  the  inlet  of  the  bellmouth  connecting  the  test  facility  plenum  to  the  engine 
inlet  must  have  the  correct  total  pressure,  total  temperature,  and  humidity  such  that  the 
condition  of  the  water-air  mixture  at  the  engine  compressor  face  simulates  the  flight 
condition.  In  general,  the  water  loading  factor  (fc)  should  be  established  equal  to  the 
water  loading  factor  of  the  free-stream  icing  conditions. 

SECTION  IV 

ANALYSIS  OF  MULTIPHASE  FLOW  WITH  APPLICATION 
TO  AN  ICING  TEST  CELL 


4.1  INTRODUCTION 

In  this  section,  an  analysis  of  a  typical  direct-connect  type  of  turbine  engine  icing 
test  cell  is  presented.  Through  a  number  of  simplifying  assumptions,  the  problem  is  reduced 
to  an  analysis  of  one-dimensional  flow  in  a  variable  area  duct.  The  analysis  provides  a 
means  for  predicting  the  thermodynamic  and  kinetic  state  of  the  air-water  mixture  at 
the  test  section.  It  is  useful  in  determining  the  test  cell  inlet  conditions  required  to  simulate 
a  given  icing  condition. 

The  analysis  parallels  the  analysis  of  exhaust  gas  spray  coolers  by  Pelton  and  Willbanks 
(Ref.  11). 

4.2  THE  BASIC  ASSUMPTIONS  USED  FOR  DEVELOPING  THE  MATHEMATICAL 

MODEL  OF  THE  MULTIPHASE  AIR-WATER  DUCT  FLOW 

Figure  8  shows  a  schematic  of  the  icing  test  facility  being  considered.  The  duct  has 
an  arbitrary  cross-sectional  area,  and  water  is  injected  parallel  to,  and  in  the  direction 
of,  the  gas  flow  at  some  location  along  the  duct.  The  analysis  is  based  on  overall 
conservation  of  species,  mass,  momentum,  and  energy  both  globally  and  in  each  phase. 
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The  mathematical  model  for  flow  in  the  duct  is  based  on  the  following  assumptions: 

1.  All  gases  including  water  vapor  are  treated  as  ideal  gases. 

2.  One-dimensional  flow. 

3.  The  gas  phase  is  homogeneous  at  any  axial  station  except  in  the  boundary 
layers  of  the  droplets. 

4.  The  drop-size  distribution  of  the  spray  can  be  characterized  by  discrete 
drop  sizes. 

5.  Drops  are  injected  parallel  to  and  in  the  direction  of  the  gas  How  and 
remain  entrained  in  the  gas  stream  throughout  the  duct.  That  is, 
de-entrainment  by  gravity  and  turbulent  diffusion  is  considered  negligible. 

6.  Each  drop  injected  maintains  its  identity  through  the  duct.  Hence,  collisions, 
agglomeration,  and  drop  breakup  are  considered  negligible. 

7.  The  drops  are  uniformly  distributed  at  any  cross  section  of  the  duct. 

8.  Internal  resistance  to  heat  transfer  within  a  drop  is  much  lower  than  external 
or  film  resistance  to  heat  transfer.  Thus,  the  internal  temperature  of  a  drop 
is  uniform  at  any  instant  of  time. 

9.  Heat  transfer  and  friction  at  the  duct  walls  are  negligible. 

10.  Cross-sectional  area  of  the  duct  is  a  prescribed  function  of  axial  location. 

11.  The  drops  are  spherical  in  shape. 

12.  Drops  can  be  supercooled  to  a  critical  temperature  which  is  a  prescribed 
function  of  drop  diameter. 

13.  Upon  reaching  the  critical  temperature,  the  drop  temperature  suddenly  rises 
to  the  freezing  point  and  becomes  a  homogeneous  equilibrium  mixture  of 
solid  and  liquid.  The  drop  remains  at  the  normal  freezing  point  until  enough 
heat  has  been  removed  to  completely  freeze  the  drop. 

14.  Vaporization  and  condensation  are  occurring  at  any  instant  at  equilibrium 
rates.  That  is,  at  any  instant,  the  drop  surface  vapor  pressure  is  equal  to 
the  saturation  pressure  corresponding  to  drop  temperature.  The  droplets 
of  interest  in  icing  problems  are  large  enough  for  the  effect  of  surface 
tension  on  vapor  pressure  to  be  negligible. 
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4.3  THE  BASIC  EQUATIONS  OF  CONSERVATION 

Based  on  the  assumptions  of  the  previous  section,  the  following  global  conservation 
equations  can  be  written: 

Conservation  of  Mass 


wv  +  Wnc  +  S  Wc.  -  W  -  Constant  (16) 

Conservation  of  Water  (Condensed  Phase  +  Vapor) 

wv  +  ?WC  -  Ww  -  Constant  (17) 

Conservation  of  Noncondensable  Gas 


®nc  **  V^l-C^A  -  Constant 


(18) 


Conservation  of  Momentum 


dC 

;2 _ I 

8  dx 


dV 


+  (l-Cv)Vg 


dx 


/  dfc  dVc\  1 

+  vf(i-c^*x(v0|-j^  +  i0i-jt)+  - 


-  c. 


ip 

dx 


=  0  (19) 


Conservation  of  Energy 


Mass,  momentum,  and  energy  transfer  between  phases  can  be  treated  by  first 
considering  these  processes  for  a  single  droplet.  For  droplet  of  class  i,  the  following 
conservation  equations  apply: 

Conservation  of  Mass 


d«d 
d  x 


k^df^ 


Conservation  of  Momentum 

dV_ 


dM, 


M 


,  V  - -  +  /v  -v)v  - - 

di  ci  dx  \  ci  8/  ci  d* 


(21) 


(22) 
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Conservation  of  Energy 


dEc  dMd 

'  ■  MflT.  -  (T.),]  +  kvlT  VCi-j-l 


c;  dx 


By  considering  all  droplets  of  class  i, 


d{c. 

fc. 

dMd 

1 

1 

1 

TT 

=  \ 

dx 

d(Kc  W  ) 

C;  Cj 

n 

ci 

dE, 

1 

dx 

Md; 

dx 

The  equation  of  state  is 


g  g 


*vMv+tl~  xy)Mnc 


The  molar  fraction  and  mass  fraction  of  vapor  are  related  by 

_Cv/Mv 

XV  =  Cv/Mv+U-.Cv)/Mnc 

Droplet  diameter  and  mass  are  related  by  the  expression: 


(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


The  molar  fraction  of  vapor  at  a  droplet  surface  is  related  to  saturation  pressure  and 
static  pressure  by 

(  v9)i  =  Pv!t,  /p  (29) 

t 

For  low  rates  of  mass  transfer,  the  heat  and  mass  transfer  coefficients  for  spherical 
droplets  can  be  obtained  from  the  semi-empirical  equation  of  Ranz  and  Marshall  (Ref. 
12).  Thus, 


13 


AEDC-TR-73-144 


k 


x 


-£=(2.0  +  0.6  Rc1/2Sc1/3) 
dScM 


(30) 


and 


F  -  -r-J-  (2.0  +  0.6  Re1/2Pr1/3) 
d  Pr 


(31) 


where  the  Reynolds  number  is  given  by 


Re 


(32) 


Carlson  and  Hoglund  (Ref.  13)  determined  the  following  semi-empirical  equation  for  the 
drag  coefficient  of  a  solid  sphere  in  low-speed  flow  and  in  the  absence  of  mass  transfer: 


CD  =  |i(l  +0.15  Re0-687) 


(33) 


This  relation  is  assumed  to  be  valid  for  low  rates  of  surface  mass  transfer  as  well.  The 
above  transport  equations  are  used  in  the  present  study  although  the  rates  of  mass  transfer 
immediately  following  droplet  injection  may  be  high.  Since  the  rate  of  mass  transfer  drops 
off  to  a  low  value  very  rapidly,  the  error  introduced  by  omitting  corrections  to  F,  kx, 
and  Cp  for  high  rates  of  mass  transfer  is  believed  to  be  small. 

4.4  A  MODEL  FOR  DROPLET  FREEZING 

Many  experiments  have  confirmed  the  fact  that  liquid  water  can  exist  in  a  supercooled 
state  (Ref.  14).  However,  the  results  of  the  many  studies  are  not  completely  in  agreement 
with  regard  to  the  parameters  governing  the  maximum  degree  of  supercooling.  Factors 
affecting  maximum  supercooling  to  some  degree  include:  rate  of  cooling,  concentration 
of  impurities,  and  drop  size  (Refs.  14  through  16). 

In  the  present  study,  the  following  simple  model  is  proposed  for  freezing  a  droplet. 
It  will  be  assumed  that  the  maximum  degree  of  supercooling  is  a  known  function  of 
only  the  droplet  diameter.  Thus,  a  drop  of  given  diameter  can  be  supercooled  to  some 
critical  temperature.  Upon  reaching  the  critical  temperature,  the  temperature  of  the  drop 
suddenly  rises  to  the  normal  freezing  point  and  becomes  a  homogeneous  equilibrium 
mixture  of  solid  and  liquid.  The  relative  fractions  of  solid  and  liquid  are  governed  by 
the  initial  degree  of  supercooling.  The  drop  remains  at  the  freezing  point  until  enough 
heat  has  been  removed  to  completely  freeze  the  drop. 
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A  curve  of  critical  temperature  as  a  function  of  drop  diameter  was  deduced  from 
data  given  in  Ref.  15  and  is  illustrated  in  Fig.  9.  The  data  can  be  correlated  by  the 
following  equations: 


Tsf  =  463.2  +  16.32  en  (0.0333d;  d  <  300  (34) 


and 

Tsf  =  463.2 


300  <  d  <  1000 


where  TSf  is  in  °R  and  d  is  the  droplet  diameter  in  microns.  This  relationship  represents 
the  average  critical  temperature  for  a  large  number  of  drops  of  given  diameter.  The  variation 
in  critical  temperatures  of  individual  drops  from  the  average  may  be  quite  large.  For 
purposes  of  the  present  study,  it  is  assumed  that  all  of  the  drops  of  a  given  diameter 
freeze  upon  reaching  the  critical  temperature  in  the  manner  given  above. 

4.5  THERMODYNAMIC  PROPERTY  DATA 

It  is  necessary  to  specify  the  thermodynamic  property  data  including  enthalpy,  specific 
heat,  and  vapor  pressure  of  each  constituent  as  a  function  of  temperature.  The  following 
equations  serve  this  purpose. 

4.5.1  Property  Data  for  Air 

For  air,  the  enthalpy  and  specific  heat  data  given  in  Ref.  17  are  satisfactorily  correlated 
by 

C  =  0.2318  +  0.1040  x  10'4T  +  0.7166  x  10‘8T2 


and 

hnc  =  0.2318T  +  0.0520  x  10'4T  +  0.2389  x  10‘8T3  -  127.064 

where  T  is  in  °R,  CPnc  is  in  Btu/lbm-°R,  and  h  is  in  Btu/lbm.  This  correlation  is  valid 
for  400  <  T  <  1700.  The  subscript  nc  is  used  since  air  is  the  noncondensable  species 
in  the  flow. 

4.5.2  Property  Data  for  Water  Vapor 

For  water  vapor,  the  enthalpy  and  specific  heat  data  of  Ref.  18  are  satisfactorily 
correlated  by 
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CD  •=  0.4304  +  0.1678  x  10‘4T  +  0.2781  x  10'7T3 

H  V 


and 

Kv  -  0.4304T  +  0.0839  x  10‘4T2  +  0.0927  x  10‘7T3  +  806.6 

where  T  is  in  °R,  CPv  is  in  Btu/lbm-°R,  hv  is  in  Btu/lbm,  and  400  <  T  <  1700. 

4.5.3  Property  Data  for  the  Condensed  Phases 

The  specific  heat  of  the  liquid  phase  was  taken  to  be  1  Btu/lbm-°R,  and  the  specific 
heat  of  the  solid  phase  was  taken  to  be  0.485  Btu/lbm-0R. 

4.5.4  Vapor  Pressure  of  the  Condensed  Phase 

The  following  equation  is  a  reasonably  good  correlation  of  the  vapor  pressure  of 
liquid  water  as  given  in  Ref.  11  over  the  temperature  range  492  <  T  <  672: 


where  T  is  in  °R,  pv  is  in  psfa,  and  hfg  is  the  latent  enthalpy  of  evaporization 
given  by 

hfg  -  1352.3  -  0.5696T  +  0.0839  x  10'4T2  +  0.0927  x  10*7T3 

The  vapor  pressure  in  psfa  of  supercooled  liquid,  that  is,  liquid  at  a  temperature 
less  than  492°R,  is  given  by  the  following  equation  as  taken  from  Ref.  14: 


Pv  = 


+ 


Cj(t2—  K) 

t 


r  (d,u2-k)2>  i  „  r 

[10  1  -  l|  +  Ejho 


(F  j  ( 37  4 . 1 1— t) 


where 

t  (=  T/1.8)  is  the  absolute  temperature  in  degrees  Kelvin,  and  Aj,  Bj,  Cj,  Dj  ,Ej,  Fj, 
and  K  are  constants  whose  values  are: 

Aj  -  5.4266514 
Bj  =  -2005.1 
Cj  =  1.3869  x  10‘4 
Dj  =  1.1965  x  10'11 
Ei  =  4.4  x  10*3 
F,  =  -5.7148  x  10'3 
Kj  =  2.937  x  10s 
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The  equation  for  the  vapor  pressure  of  ice  in  psfa  as  given  in  Ref.  14  is 
pv  =  2.7845  exp  [2.3[A2/t  +  0.4343B2  Mt)  +  C2t  +  D2t2  +  E2]l 


where 


t  (=T/1.8)  is  the  absolute  temperature  in  degrees  Kelvin  and  A2,  B2,  C2,  D2, 
and  E2  are  constants  whose  values  are: 

A 2  =  -2.4455646  x  103 
B2  =  8.2312 
C2  =  1.677006  x  10'2 
D2  -  1.20514  x  10'5 
E2  =  -6.757169 

4.6  THE  COMPUTER  PROGRAM 

The  system  of  equations  is  closed  and  can  be  integrated  for  specified  initial  conditions 
to  give  the  variation  in  thermodynamic  and  kinetic  states  along  the  duct.  The  equations 
were  programmed  in  FORTRAN  IV  language  for  solution  on  a  digital  computer.  A  program 
listing  is  given  in  Appendix  III.  The  computer  program  developed  in  connection  with  a 
study  of  spray  coolers  and  reported  in  Ref.  1  1  was  used  as  the  basis  for  the  program 
in  Appendix  III. 

The  program  has  the  capability  to  handle  multiple  water  spray  stations  at  arbitrary 
locations  along  the  duct.  The  drops  originating  at  any  station  must  have  uniform  diameters. 
The  drop  diameter  may,  however,  vary  from  station  to  station.  With  this  capability,  a  spray 
with  a  distribution  of  drop  sizes  can  be  treated,  providing  the  distribution  can  be 
characterized  by  a  finite  number  of  discrete  drop  sizes.  To  accomplish  this,  each  drop-size 
class  (i)  is  assigned  to  an  individual  injection  station  (i).  From  a  given  distribution  of 
drop  sizes,  the  fraction  of  injected  liquid  containing  droplets  in  size  class  (i)  can  be 
determined.  This  fraction  corresponds  to  the  fraction  of  the  total  liquid  injection  rate 
to  be  assigned  to  injection  station  (i). 

In  the  computer  program,  droplet  freezing  is  handled  in  the  following  way:  The 
temperatures  of  drops  originating  at  each  injection  station  are  monitored,  and  when  the 
temperature  of  a  droplet  reaches  its  freezing  temperature  defined  by  Eq.  (34),  then  the 
droplet  is  considered  to  start  freezing.  An  energy  balance  predicts  the  state  of  the  droplet 
at  each  step  in  the  flow  while  the  droplet  temperature  is  held  constant  at  492°R  (32°F). 
After  sufficient  energy  is  lost  from  the  droplet  by  conduction  and  evaporation  to  turn 
the  droplet  completely  to  ice,  the  droplet  model  is  replaced  by  an  ice  particle  model, 
and  its  temperature  is  allowed  to  vary.  For  the  rest  of  the  calculations,  the  particular 
size  drop  that  has  frozen  is  treated  as  a  spherical  particle  of  ice,  and  its  temperature 
is  again  determined  by  energy  and  mass  transfer  at  the  ice  particle  surface. 
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4.7  SUMMARY 

An  analysis  of  the  multiphase  flow  phenomena  occurring  in  a  typical  turbine  engine 
icing  test  cell  duct  flow  was  presented.  The  analysis  was  based  on  one-dimensional  flow 
of  a  mixture  of  air,  water  vapor,  liquid  water,  and  ice  in  a  variable  area  duct. 

In  the  following  section,  the  analysis  and  computer  program  are  used  to  make 
calculations  for  a  typical  icing  test  cell.  A  parametric  study  is  performed  to  determine 
the  influence  of  the  initial  conditions  on  the  final  thermodynamic  and  kinetic  states. 

SECTION  V 

APPLICATION  OF  THE  MULTIPHASE  FLOW 
ANALYSIS:  A  PARAMETRIC  STUDY 


5.1  INTRODUCTION 

In  this  section,  results  obtained  with  the  multiphase  flow  analysis  for  flow  in  an 
icing  test  facility  under  typical  operating  conditions  are  reported  and  discussed.  A 
parametric  study  was  performed  to  determine  the  influence  of  the  inlet  conditions  on 
the  kinetic  and  thermodynamic  flow  properties  at  the  test  section. 

The  particular  duct  geometry  chosen  for  the  study  is  shown  in  Fig.  10.  The  bellmouth 
entrance  approximates  the  inlet  section  of  the  AEDC  Propulsion  Engine  Test  Cel]  (J-I ) 
and  typifies  many  other  installations.  The  inlet  parameters  which  were  used  are: 

1.  The  primary  thermodynamic  and  kinetic  state  variables  such  as  pressure, 
temperature,  bellmouth  inlet  velocity,  and  humidity 

2.  The  water  loading  factor  (fc)  which  gives  the  mass  of  liquid  water  sprayed 
into  the  airflow  at  the  duct  inlet  per  unit  mass  of  air 

3.  The  water  injection  velocity 

4.  The  mass  median  water  droplet  diameter 

5.2  THE  EFFECT  OF  MATHEMATICALLY  SIMULATING  THE  WATER  SPRAY 

CLOUD  WITH  A  DISTRIBUTION  OF  DROPLETS  AS  COMPARED  WITH 

SIMULATING  THE  WATER  SPRAY  CLOUD  WITH  A  SINGLE-SIZE 

DROPLET 

The  analysis  was  first  applied  to  study  the  effect  of  a  spray  cloud  with  a  droplet 
distribution  on  the  kinetic  and  thermodynamic  state  of  flow  through  the  duct  in 
comparison  with  the  effects  of  a  spray  having  only  a  single-sized  drop.  For  the  latter 
spray,  the  drop  size  was  equal  to  the  mass  median  diameter  of  the  given  distribution. 
The  distribution  was  obtained  from  Ref.  9  and  is  a  measured  distribution  from  an  icing 
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simulation  test  utilizing  such  a  duct  flow.  The  mass  median  diameter  was  approximately 
25  microns,  and  the  mass  mean  diameter  was  approximately  18  microns.  Figure  11  shows 
three  drop-size  distributions  that  were  presented  in  Ref.  9.  The  distribution  utilized  in 
this  study  is  distribution  A-3  of  Fig.  11.  The  specific  purpose  of  this  mathematical 
experiment  was  to  determine  theoretically  if  the  practice  of  characterizing  an  engine  icing 
condition  by  the  liquid  water  content  of  the  flow  and  a  mean  effective  drop  diameter 
in  the  spray  cloud  is  a  valid  method  for  correlating  the  final  state  of  the  water-airflow 
at  the  engine  station  with  the  given  duct  inlet  conditions.  If  not,  then  it  is  necessary 
to  take  into  account  the  actual  distribution  of  droplet  sizes  produced  in  the  spray  cloud 
by  the  atomizing  nozzle  system  when  trying  to  correlate  engine  icing  conditions  with 
duct  inlet  conditions.  Thus,  another  possible  effect  of  the  droplet  distribution,  other  than 
on  the  capture  efficiency  of  the  component  parts  of  the  engine  inlet,  may  be  to  affect 
the  thermodynamic  and  kinetic  state  of  the  multiphase  duct  flow. 

The  results  of  the  calculations  made  with  the  droplet  distribution  mathematical 
simulation  and  the  mean  effective  diameter  simulation  are  shown  and  compared  in  Figs. 
12,  13  and  14.  The  other  initial  conditions  assumed  are  presented  in  these  figures.  In 
Fig.  12,  the  local  static  gas  temperature  (Tg)  is  plotted  versus  length  from  the  duct  inlet. 
The  gas  temperature  is  not  highly  dependent  on  whether  the  spray  cloud  is  characterized 
by  having  a  distribution  of  droplet  sizes  or  a  single  uniform  drop  size.  Figure  13  shows 
the  liquid  droplet  temperatures  for  each  size  drop  in  the  distribution  plotted  against  length 
from  the  duct  inlet.  Again,  while  the  larger  droplets  in  the  spray  cloud  required  greater 
length  of  flow  before  reaching  a  quasi-equilibrium  temperature,  all  size  droplets  have 
reached  approximately  the  same  temperature  at  the  engine  station,  and  this  is  the  same 
temperature  reached  by  the  mass  median  drop  (Fig.  14).  Therefore,  for  typical  icing 
conditions,  the  characterization  of  the  icing  conditions  by  liquid  water  content  and  mean 
effective  drop  diameter  is  sufficiently  accurate  for  testing  purposes.  The  engine  icing 
conditions  can  be  correlated  in  terms  of  these  two  parameters. 

By  using  the  droplet  freezing  model  in  the  program,  a  comparison  was  made  with 
the  same  model  droplet  distribution  spray  cloud  mentioned  earlier  and  the  single  uniform 
size  spray  cloud  of  the  mean  effective  diameter  of  the  distribution.  The  computations 
showed  that  for  the  typical  icing  test  conditions,  the  freezing  model  chosen,  and  droplet 
sizes  for  icing,  there  was  no  droplet  freeze-out.  The  droplets  remain  subcooled  and  unfrozen 
for  duct  lengths  of  up  to  40  ft  from  inlet  to  engine  station.  Thus,  the  icing  conditions 
defined  by  the  military  specifications  for  engine  testing  can  be  characterized  by  liquid 
water  content  and  mass  median  diameter,  and  for  these  conditions  there  should  be  no 
change  in  the  liquid  water  content  because  of  droplet  freeze-out.  It  should  be  noted, 
of  course,  that  in  a  real  test  facility,  droplet  freeze-out  may  occur  when  there  is 
contamination  of  the  spray  water  with  particles  that  enable  ice  nucleation  and  crystal 
formation.  However,  this  problem  was  not  considered  in  this  study. 
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5.3  PARAMETRIC  STUDY  OF  THE  EFFECTS  OF  THE  THERMODYNAMIC 

AND  KINETIC  DUCT  INLET  PARAMETERS  ON  THE  STATE  OF  THE 

FLOW  AT  THE  TEST  SECTION 

5.3.1  Description  of  the  Baseline  Values  of  the  Duct  Inlet  Parameters 

For  the  duct  geometry  previously  discussed,  a  parametric  study  was  carried  out  by 
varying  the  duct  inlet  properties  one  at  a  time.  The  inlet  variables  that  were  studied  are 
given  below  along  with  their  baseline  values.  The  baseline  values  of  the  duct  inlet  (initial) 
parameters  are  as  follows: 

1.  Liquid  water  initial  temperature,  Ts  (632°R  or  172<>F) 

2.  Liquid  water  content,  L  (1  gm/m3  or  6.234  x  10* 5  lbm/ft3) 

3.  Inlet  air  specific  humidity,  (3.715  x  10*3  lbm  f^O/lbm  air) 

4.  The  mass  median  drop  diameter,  dm  (3.175  x  1 0~4  ft  or  96.8  microns) 

5.  The  water  injection  velocity,  Vs  (8.35  ft/sec) 

6.  Air  inlet  velocity,  Vg  (80.6  ft/sec) 

7.  Duct  radius  at  constant  area  section,  R  (3.8  ft) 

The  reason  dm  was  chosen  at  approximately  100  microns  rather  than  a  more 
representative  value  of,  say,  25  microns  is  simply  that  small  drops  cool  faster  than  larger 
drops.  Thus,  for  the  larger  baseline  droplet  chosen,  the  axial  variations  of  the  other 
properties  were  more  clearly  defined  when  the  baseline  drop  diameter  was  100  microns 
than  if  it  were  25  microns.  This  effect  will  be  evident  in  the  figure  to  be  presented  showing 
effects  of  varying  drop  diameter.  The  results  of  the  parametric  study  are  discussed  in 
the  following  sections. 

5.3.2  The  Effect  of  Variation  of  Initial  Water  Temperature 

The  effect  of  the  temperature  variation  of  the  water  spray  into  the  bellmouth  was 
negligible;  that  is,  the  final  state  of  the  air-water  mixture  at  the  engine  station  did  not 
depend  significantly  on  the  temperature  of  the  injected  water.  These  results  are  shown 
in  Figs.  15  and  16  where  the  temperature  difference  (AT)  is  plotted  versus  length  from 
the  duct  entrance.  The  obvious  explanation  of  these  results  is  that  the  amount  of  water 
being  sprayed  into  the  bellmouth  in  an  icing  simulation  (and  hence  the  amount  of  energy 
carried  into  the  bellmouth  by  the  water)  is  extremely  small  compared  with  the  mass  flow 
(or  the  energy  flow)  of  air.  Thus,  the  process  in  the  duct  for  icing  simulation  is  nearly 
a  constant  energy  flow,  not  considering  heat  transferred  through  the  duct  walls.  The  gas 
temperature  at  the  engine  station  depends  primarily  on  the  gas  total  temperature  and 


20 


AEDC-TR-73144 


the  Mach  number  of  the  flow  at  the  engine  station.  The  water  droplets  cool  quickly  down 
to  approximately  the  same  quasi-equilibrium  temperature  which,  it  was  found  later,  could 
be  predicted  in  terms  of  the  air  temperature  and  humidity  at  the  engine  station.  It  should 
be  emphasized  that  this  conclusion  applies  only  for  small  liquid  water  loading  typical 
of  icing  simulations. 

5.3.3  The  Effect  of  Variation  of  Initial  Liquid  Water  Content 
(Loading  Factor) 

The  effect  of  the  variation  of  the  liquid  water  content  (L)  from  about  1  to  10  gm/m3 
on  the  Final  temperature  of  the  airflow  is  shown  in  Fig.  17  to  be  small.  Again,  the  reason 
for  this  small  effect  is  that  the  absolute  quantity  of  water  sprayed  into  the  duct  flow 
is  very  small  compared  with  the  mass  flow  of  air.  The  variation  of  L  from  1  to  10  gm/m3 
caused  only  a  minor  variation  in  the  gas  temperature  at  the  engine  station  (4°R)  and 
extremely  small  effect  on  pressure  and  velocity  of  the  airflow  at  the  engine  station.  Also, 
the  temperature  of  the  liquid  water  at  the  engine  station  did  not  change  significantly 
with  variation  of  the  liquid  water  content.  This  is  shown  in  Fig.  18  where  AT  is  plotted 
versus  duct  length.  In  this  figure,  there  is,  at  most,  two  degrees  difference  between  the 
water  temperature  at  the  engine  station  over  this  range  of  liquid  water  content. 

5.3.4  The  Effect  of  Variation  in  Mass  Median  Droplet  Diameter 

The  effect  of  varying  the  mean  effective  diameter  is  shown  in  Fig.  19.  As  stated 
previously  and  shown  in  Fig.  19,  small  drops  cool  down  faster  than  large  drops.  As  Fig. 
19  indicates,  a  spray  cloud  of  uniform  sized  drops  of  25-micron  diameter  cools  down 
at  a  much  faster  rate  than  a  spray  cloud  of  100-micron  diameter  drops.  The  25-micron 
drop  temperature  drops  below  the  gas  temperature  but  then  recovers  and  approaches  the 
gas  temperature.  Droplets  of  200  microns  and  larger  will  be  substantially  out  of  equilibrium 
with  the  gas  temperature  at  25  ft  from  the  duct  inlet.  This  size  of  droplets  requires  a 
duct  length  substantially  longer  than  25  ft  before  their  temperature  would  come  to  within 
a  few  degrees  of  the  gas  temperature.  The  water  spray  system  must  produce  only  a  very 
small  number  of  drops  larger  than,  say,  200  microns  to  minimize  their  effect  on  the  icing 
simulation. 

5.3.5  The  Effect  of  Variation  in  Inlet  Air  Humidity 

The  effects  of  varying  the  inlet  air  relative  humidity  from  zero  (completely  dry  air) 
to  100  percent  (completely  saturated  air)  are  shown  in  Fig.  20.  It  can  be  noted  that 
the  inlet  humidity  is  very  important  in  determining  the  temperature  of  the  water  droplets. 
The  reason  is  that  evaporation  and  condensation  at  the  droplet  surface  play  major  parts 
in  determining  droplet  temperature.  In  general,  the  rate  of  surface  mass  transfer  is  highly 
dependent  on  the  humidity  of  the  air  surrounding  the  droplet.  By  controlling  the  inlet 
humidity,  it  is,  therefore,  possible  to  have  some  degree  of  control  over  the  thermal  state 
in  the  test  section. 
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5.3.6  The  Effect  of  the  Difference  between  Inlet  Airflow  Velocity 
and  Inlet  Water  Velocity 

The  effect  of  the  variation  of  the  water  injection  velocity  on  the  water  temperature 
at  the  engine  station  is  shown  in  Fig.  21.  This  Figure  shows  that  as  the  velocity  difference 
between  the  water  spray  cloud  and  the  airflow  is  decreased,  a  greater  length  of  flow  is 
required  before  the  water  temperature  approaches  the  air  temperature;  that  is,  increasing 
the  relative  velocity  increases  the  heat  and  mass  transfer  from  the  drops.  The  coefficients 
of  mass  and  heat  transfer  are  defined  by  Eqs.  (24)  and  (25),  respectively.  These  coefficients 
are  dependent  on  a  Reynolds  number  based  on  the  relative  velocity  or  velocity  difference 
between  the  airflow  and  the  water  drops.  Thus,  heat  and  mass  transfer  increase  as  the 
relative  velocity  between  air  and  the  water  drops  increases.  Hence,  control  of  the  initial 
velocity  difference  between  the  sprayed  water  and  the  airflow  at  the  duct  inlet  may  provide 
a  method  for  controlling  the  flow  conditions  at  the  engine  face.  Note  that  this  comparison 
was  computed  with  dm  =  15  microns  because  the  effect  of  inlet  velocity  difference  was 
more  pronounced  than  a  calculation  made  with  dm  =  96.8  microns. 

5.3.7  Discussion 


After  injection  of  a  droplet  into  an  air  stream,  the  velocity  of  the  droplet  tends 
to  approach  the  velocity  of  the  air.  The  temperature  of  the  droplet  approaches  the  wet 
bulb  temperature  of  the  air;  hence,  unless  the  air  is  saturated,  the  droplet  and  air 
temperature  will  not  become  equal.  In  a  duct  flow,  equilibrium'  can  be  achieved  in  a 
sufficiently  long  duct  since  droplet  evaporation  or  condensation  leads  to  saturation  of 
the  air.  The  duct  length  for  this  to  occur  is  several  orders  of  magnitude  greater  than 
the  length  available  for  a  practical  test  cell.  Figures  12  through  17  show  this  effect.  It 
can  be  noted  that  the  droplet  temperature  usually  becomes  essentially  constant  only  a 
few  feet  downstream  of  the  injection  station.  This  distance  is,  of  course,  a  function  of 
the  inlet  flow  conditions  and  the  droplet  diameter. 


The  steady-state  droplet  temperature  (Ts)  (wet  bulb  temperature)  can  be  calculated 
from  the  relation: 


T  -  T  =  -hf 

g  S  I 

r  P 


P r  Mv  /  |  /  \ 

Sc  MC  \  v  Pv  T,  V 


where  Tg  is  the  local  gas  temperature,  Cp  is  the  specific  heat  of  air,  xv  is  the  mole 
fraction  of  water  vapor  in  the  gas  phase,  PvIts  *s  the  saturation  vapor  pressure  of  water 
at  temperature  Ts,  and  p  is  the  local  static  pressure.  This  equation  can  be  used  to  estimate 
the  temperature  nonequilibrium  in  the  test  section.  An  accurate  value  of  xv  can  be  obtained 
only  by  calculating  the  entire  duct  flow.  However,  for  many  flow  conditions  of  interest 
in  icing  simulation,  the  specific  humidity  (cj),  and  therefore  xv,  is  essentially  constant 
along  the  duct.  Thus  for  an  estimate,  xv  can  be  taken  to  be  equal  to  the  vapor  content 
upstream  of  the  spray  bank. 
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It  is  very  difficult  to  obtain  an  accurate  measurement  of  the  humidity  under  icing 
conditions.  If  accurate  knowledge  and  control  of  humidity  is  important  in  an  icing  test, 
it  may  be  desirable  to  circumvent  the  difficult  problem  of  humidity  measurement  by 
creating  known  humidity  in  the  test  cell.  This  can  be  accomplished  by  supplying  the  test 
cell  with  very  dry  air  and  mixing  the  required  amount  of  superheated  steam  with  the 
air  to  create  the  desired  humidity.  It  is  important  for  the  mixing  to  be  essentially  complete 
by  the  time  the  flow  reaches  the  water  sprays.  Very  good  humidity  control  should  be 
possible  with  this  system. 

5.3.8  Summary 

A  summary  of  the  results  of  this  section  of  the  report  follows.  First,  for  the  flow 
conditions  assumed,  it  was  found  that  the  thermodynamic  and  kinetic  state  of  the  air 
and  water  mixture  at  .  the  test  section  of  a  typical  direct-connect  icing  facility  duct  is 
not  dependent  on  the  drop-size  distribution.  The  analysis  showed  that  after  about  25 
ft  downstream  of  the  inlet,  the  spray  cloud  can  be  satisfactorily  modeled  with  two 
parameters,  the  mass  median  drop  diameter  and  the  liquid  water  content.  The  droplet 
size  distribution  used  in  the  calculation  was  a  measured  distribution  from  an  icing  test 
cell.  Second,  the  temperature  of  the  air,  water  vapor,  and  water  mixture  25  ft  downstream 
of  the  inlet  was  practically  independent  of  the  initial  temperature  of  the  spray  water. 
The  reason  is  that  the  liquid  water  contents  are  extremely  small.  Third,  the  effect  of 
varying  the  liquid  water  content  between  I  and  10  on  the  thermal  and  kinetic  states 
at  the  test  section  is  small.  Fourth,  the  effect  of  inlet  air  humidity  upstream  of  the  water 
spray  station  is  relatively  important  in  determining  the  test  section's  thermal  state  and 
a  substantial  degree  of  control  over  the  thermal  state  at  the  test  section  is  possible  through 
humidity  control.  Fifth,  the  temperature  of  relatively  small  (5  to  25  microns)  droplets 
approaches  the  air  temperature  at  a  greater  rate  than  that  of  the  relatively  large  (100 
to  1000  micron)  droplets  because  their  mass  is  considerably  less.  Hence,  the  small  droplets 
are  effectively  at  the  air  temperature  in  a  much  shorter  length  of  duct  than  the  relatively 
large  droplets,  and  for  this  reason,  the  test  facility  water  spray  systems  must  not  produce 
a  significant  fraction  of  water  drops  greater  than  100  microns  if  effective  simulation  is 
to  be  provided.  Sixth,  the  water  spray  velocity  should  be  substantially  different  from 
the  surrounding  airflow  velocity  in  order  to  increase  the  rate  of  heat  and  mass  transfer 
from  the  drop  and  the  approach  to  equilibrium.  That  is,  a  large  initial  degree  of  kinetic 
nonequilibrium  between  the  injected  liquid  and  air  helps  to  promote  a  more  rapid  approach 
to  thermal  equilibrium. 


SECTION  VI 

CONSIDERATION  OF  SOME  FACTORS  AFFECTING 
ICING  SIMULATION 


6.1  INTRODUCTION 

Since  most  ground  test  facilities  cannot  duplicate  all  of  the  flow  conditions 
experienced  by  an  engine  in  flight,  it  is  necessary  to  determine  which  of  the  icing  variables 
are  relatively  most  important  in  determining  engine  icing  behavior.  In  this  section  of  the 
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report,  the  effect  of  droplet-size  distribution  on  capture  efficiency  is  discussed.  In  addition, 
an  analysis  of  the  heating  rate  required  to  maintain  an  anti-icing  condition  on  a  surface 
in  an  icing  environment  is  provided.  The  relative  importance  of  correctly  simulating  the 
specific  humidity,  liquid  water  content,  and  mass  median  droplet  diameter  for  simulating 
icing  conditions  in  a  propulsion  engine  test  facility  is  illustrated  quantitatively. 

6.2  THE  EFFECT  OF  DROPLET-SIZE  DISTRIBUTION 

ON  CAPTURE  EFFICIENCY 

In  Section  V,  the  effects  of  a  water  spray  cloud  with  a  distribution  of  droplet  sizes 
were  compared  with  the  effects  of  a  spray  cloud  with  a  uniform  droplet  size  on  the 
final  thermodynamic  and  kinetic  state  of  a  flowing  mixture  of  water  and  air.  The  effects 
were  shown  to  be  equivalent  for  a  typical  icing  situation. 

There  are,  however,  other  possible  effects  which  the  drop-size  distribution  might  have 
which  must  be  considered  when  simulating  icing  in  ground  test  facilities.  For  example, 
a  given  drop-size  distribution  might  affect  ice  accumulation  rates  on  the  exposed  surfaces 
of  components  of  the  engine  inlet  by  affecting  the  capture  efficiency  of  these  surfaces. 
Capture  efficiency  is  defined  here  as  follows.  Consider  the  volume  swept  by  the  projected 
area  of  a  surface  moving  through  a  uniform  stream.  For  a  uniform  stream  laden  with 
particles,  the  hypothetical  maximum  particle  impingement  rate  is  simply  the  product  of 
the  particle  concentration  per  unit  volume  times  the  volume  flow  swept  by  the  surface; 
that  is, 


w 

"IMF,  max 


x 


(35) 


where  Ap  is  the  projected  area  of  the  surface.  In  general,  however,  since  the  surface  disturbs 
the  flow  field  ahead  of  it,  the  particles  lying  in  the  path  of  the  surface  will  be  influenced 
by  the  disturbed  flow.  Depending  on  particle  inertia,  the  particles  will  follow  or  at  least 
be  sufficiently  influenced  by  the  streamlines  of  the  airflow  to  take  curved  trajectories, 
and  some  particles  in  the  path  of  the  surface  will  pass  over  the  surface  without  impinging. 
The  two  possibilities  are  shown  in  Fig.  22;  the  larger  droplets  in  the  stream  will  be  least 
likely  to  follow  streamlines  because  of  their  inertia.  It  is  possible  to  define  the  absolute 
rate  of  impingement  of  particles  on  a  surface  by  analyzing  the  action  of  the  particles 
in  the  disturbed  flow  field  about  the  surface.  For  a  given  size  particle,  a  limiting  trajectory 
can  be  found  which  is  just  tangent  to  the  surface.  The  tangential  trajectory  is  used  to 
find  the  cross-sectional  area  of  the  streamtube  far  upstream  of  the  surface  which  contains 
all  of  the  particles  of  a  given  size  that  will  impinge  on  the  surface.  Particles  outside  this 
streamtube  will  be  swept  over  the  surface  and  will  not  impinge.  Hence,  the  actual 
impingement  rate  is 


ftlMP  _  L»o  x  AdT  x  Voo  (36) 

where  AdT  is  the- cross-sectional  area  of  the  streamtube  defined  by  the  tangential  trajectory 
(Fig.  22).  This  analysis  assumes  that,  for  the  small  stream  loadings  typical  of  icing 
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conditions,  the  flow  field  over  the  surface  is  not  influenced  or  modified  by  the  presence 
of  the  water  droplets.  The  flow  field  is  determined  only  by  the  surface  geometry. 

The  capture  efficiency  of  a  surface  can  now  be  defined  as 

r  S'MP 

Ec  =  ft 

IM  P,  max 

For  a  given  surface  in  an  icing  condition,  each  size  water  droplet  in  the  free  stream  will 
have  a  characteristic  tangential  trajectory  because  of  its  inertia;  the  capture  efficiency  (Ec) 
of  a  surface  refers  to  a  specified  drop  size.  The  overall  capture  efficiency  of  a  body  in 
a  stream  laden  with  particles  depends  on  the  distribution  of  drop  sizes  in  the  stream. 
It  is  of  interest  to  this  study  to  determine,  at  least  qualitatively,  if  the  overall  capture 
efficiencies  of  the  engine  inlet  components  are  sensitive  to  drop-size  distributions.  If  so, 
then  this  may  be  an  important  parameter  to  simulate,  or  duplicate  if  possible,  in  the 
ground  test  facility  for  the  following  reason. 

If  the  capture  efficiencies  of  surfaces  within  an  engine  are  dependent  on  the  drop-size 
distribution,  then  it  is  possible  that  power  requirements  for  keeping  a  surface  inside  the 
engine  inlet  ice  free  would  depend  not  only  on  cloud  liquid  water  content  and  mass  median 
droplet  diameter,  but  also  on  the  icing  cloud  drop-size  distribution.  In  this  event,  it  would 
be  necessary  to  duplicate  both  the  mass  median  diameter  and  the  drop-size  distribution 
to  simulate  the  icing  environments  in  the  propulsion  test  facility.  For  this  reason,  the 
variation  of  capture  efficiency  of  typically  sized  surfaces  to  different  measured  drop 
distributions  was  investigated.  In  addition,  the  sensitivity  of  anti-icing  power  requirements 
of  exposed  surfaces  to  variations  in  stream  liquid  water  content  was  analyzed  assuming 
a  surface  capture  efficiency  of  unity. 

6.3  SURFACE  CAPTURE  EFFICIENCY  FOR  DROP  DISTRIBUTIONS  TYPICAL 

OF  ICING  CONDITIONS 

6.3.1  Introduction 

The  problem  is  to  estimate  capture  efficiencies  for  typically  sized  engine  inlet 
components  for  droplet  sizes  representative  of  icing  environments  and,  from  these 
calculations,  to  determine  if  these  capture  efficiencies  are  sensitive  to  droplet  distributions 
typical  of  icing  conditions  in  the  atmosphere  and  in  a  ground  test  facility.  The  principal 
inlet  components  of  the  high-bypass  turbofan-type  engine  considered  in  this  study  are 
the  spinner  bullet  nose  on  the  compressor  axis,  the  compressor  inlet  guide  vanes,  and 
the  first-stage  blades.  The  capture  efficiencies  of  the  guide  vanes  and  the  compressor  blades 
may  be  very  crudely  obtained  by  analogy  to  equivalent  cylinders.  The  thickness  of  the 
guide  vanes  and  compressor  blades  define  the  diameter  of  equivalent  cylinders.  The  spinner 
hub  capture  efficiency  can  be  estimated  directly  from  known  analytical  and  experimental 
results  for  bodies  of  revolution. 
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The  work  of  Langmuir  (Ref.  19)  can  be  used  to  determine  the  capture  efficiencies 
of  right  circular  cylinders  in  incompressible  flow,  given  the  cylinder  diameter,  the 
temperature  and  density  of  the  air,  and  the  density  of  the  droplets  and  their  diameter. 
The  effect  of  compressibility  has  been  found  to  be  negligible  (Ref.  20);  hence,  Langmuir's 
work  also  applies  with  sufficient  accuracy  to  compressible  flow  over  cylinders.  Langmuir 
defined  two  correlation  parameters  which  enable  the  capture  efficiency  of  a  given 
cylinder/airflow/drop  combination  to  be  estimated  from  his  analytical  results  presented 
in  Ref.  19.  These  two  parameters  are  the  inertia  parameter  (K)  where 


and  the  parameter  (0)  called  the  altitude  parameter  in  Ref.  19  as  defined  by 


4>  = 


9p2 D  V 
HPc 


By  using  the  given  information  on  a  particular  cylinder/airflow/drop  combination,  K  and 
4>  can  be  evaluated,  and  the  capture  efficiency  of  the  cylinder  under  the  given  conditions 
can  be  obtained  from  Fig.  1  or  2  of  Ref.  19.  Thus,  for  a  given  distribution  of  droplet 
sizes,  a  mass  fraction  (of  the  total  water  content,  per  unit  volume  of  gas)  can  be  calculated 
for  each  droplet-size  class.  From  this,  the  mass  fraction  of  liquid  water  impinging  on  the 
cylinder  from  each  droplet-size  class  can  be  determined.  Thus,  an  overall  mass-weighted 
capture  efficiency  of  a  cylinder  for  a  given  drop-size  distribution  can  be  calculated.  It 
should  be  mentioned  here  that  Langmuir  found  that,  for  his  assumed  droplet-size 
distributions,  one  can  equate  the  overall  capture  efficiency  of  a  cylinder  for  a  distribution 
of  droplets  to  the  capture  efficiency  of  the  cylinder  for  the  mass  median  droplet.  On 
the  other  hand  Dorsch,  et  al.  (Ref.  21),  show  that  the  capture  efficiency  of  a  5  to 
1  ellipsoid  of  revolution  for  a  given  drop  distribution  was  not  equal  to  the  capture  efficiency 
of  the  ellipsoid  for  the  mass  median  drop.  Thus,  to  clarify  the  problem  for  simulated 
icing  testing,  these  effects  were  reinvestigated  during  this  study. 


6.3.2  Droplet -Size  Distributions  in  Icing  Conditions 


The  first  problem  was  to  find  reported  drop-size  distributions  in  the  natural 
atmospheric  environment.  Distributions  typical  of  conditions  in  cumulus  clouds  are 
reported  by  Battan  and  Reitan  (Ref.  22)  and  in  fog  conditions  by  Eldridge  (Ref.  23). 
In  these  references,  the  distributions  were  presented  in  terms  of  frequency  of  occurrence 
(number  of  drops  of  a  given  size  per  sample)  versus  drop  diameter.  The  resulting  mass 
fraction  distributions  are  shown  in  Fig.  23.  Distribution  N-l  is  the  mass  fraction  distribution 
for  icing  conditions  in  fogs,  reported  in  Ref.  23.  Distribution  N-2  corresponds  to  conditions 
in  tropical  cumuliform  clouds,  and  distribution  N-3  corresponds  to  the  conditions  in 
cumuliform  clouds  over  the  L’SA;  both  of  these  distributions  were  reported  in  Ref.  22. 
In  this  study,  these  distributions  were  reduced  to  and  presented  in  the  form  of  mass 
fraction  of  total  liquid  water  content  (per  drop  size  class  per  sample)  versus  drop  diameter. 
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Drop  distributions  that  were  measured  in  an  icing  test  in  a  ground  test  facility  are 
reported  by  Gall  and  Floyd  (Ref.  9).  These  distributions  were  also  reduced  in  this  study 
to  mass  fraction  (of  total  liquid  water  content  per  drop-size  class)  versus  drop  diameter. 
These  distributions  are  presented  in  Fig.  1 1  of  this  report.  The  mass  median  and  mass 
mean  diameters  of  these  distributions  are  also  presented  in  Fig.  1  1  with  sufficient 
supplementary  data  to  identify  each  distribution  in  its  original  form  as  given  in  F'ig.  23 
of  Ref.  10. 

A  comparison  of  Figs.  10  and  23  of  this  report  shows  that  the  drop  distributions 
in  the  icing  facility  are  fairly  similar  to  drop  distributions  in  cumuliform  clouds.  To  compare 
these  distributions  with  the  distributions  utilized  by  Langmuir  and  Blodgett  (Ref.  19), 
the  original  forms  of  the  distributions  (number  frequency  versus  drop  diameter)  were 
replotted  in  this  study  in  the  form  of  number  frequency  of  occurrence  versus  a  reduced 
droplet  diameter  which  was  the  form  favored  by  Langmuir  and  Blodgett.  The  reduced 
diameter  (dfO  is  defined  as 


Comparisons  are  shown  in  Figs.  24  and  25.  In  Fig.  24,  the  cumuliform  distributions  of 
Ref.  22  and  the  artificial  distribution  of  Ref.  10  are  compared,  whereas  in  Fig.  25  the 
distributions  of  Ref.  22  are  compared  with  the  Langmuir  distributions  which  are  denoted 
as  B,  C,  D,  and  E  in  Ref.  19.  Because  the  Langmuir  distributions  were  suggested  by 
data  from  cloud  measurements  made  by  Houghton  and  Radford  (Ref.  24)  and  by  Radford 
(Ref.  25),  the  favorable  comparison  in  Fig.  25  is  not  surprising.  The  important  point 
to  be  made  from  these  comparisons  is  that  the  test  facility  spray  system,  like  that  described 
by  Gall  and  Floyd,  is  capable  of  simulating  the  droplet  size  and  water  mass  distributions 
which  occur  in  natural  icing  and  cloud  conditions  in  the  atmosphere. 

6.3.3  Effect  of  Drop  Distributions  on  Capture  Efficiency 

The  question  previously  posed  was  whether  the  capture  efficiency  of  typical  engine 
inlet  surfaces  was  dependent  on  drop  distributions  in  icing  conditions.  Since  natural  and 
artificial  icing  cloud  drop  distributions  were  shown  to  compare  favorably,  it  may  seem 
irrelevant  to  pursue  the  question  of  the  relation  between  capture  efficiency  and  drop 
distribution.  Nevertheless,  the  values  of  these  estimated  capture  efficiencies  are  important 
for  estimating  water  impingement  rates,  and  in  addition,  while  the  drop-size  distributions 
compared  favorably,  they  were  not  identical.  For  these  reasons,  overall  capture  efficiencies 
for  two  different  diameter  cylinders  and  a  5:1  ellipsoid  of  revolution  were  computed  for 
distributions  N-2,  N-3,  A-l,  A-2,  and  A-3. 

The  stream  conditions  were  the  same  for  each  case,  defined  by 

VM  =  440  ft/sec 

p  =  1932  psfa 
T  =  480°R 
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For  the  circular  cylinders,  the  inertia  parameter  (K)  and  the  altitude  parameter  (</>)  were 
used  together  with  the  results  in  Figs.  I  and  2  of  Ref.  19,  to  obtain,  for  each  droplet 
distribution,  a  mass-weighted  overall  capture  efficiency  (Eco).  The  results  of  these 
computations  are  presented  in  the  following  table: 

CYLINDER  DIAMETER:  0.04166  ft  (1/2  in.) 


Capture  Efficiency, 

Eco 

Distribution 

dm  > 

microns 

Based 

on  Distribution 

Based  on  Mass 
Median  Diameter 

N-2 

40.0 

0.928 

0.962 

N-3 

26.5 

0.923 

0.925 

A-l 

38.0 

0.950 

0.958 

A-2 

30.0 

0.938 

0.942 

A- 3 

22.0 

0.909 

0.912 

CYLINDER  DIAMETER 

i:  0.25  ft 

(3  in.) 

Capture  Efficiency, 

Eco 

Distribution 

dm  < 

microns 

Based 

on  Distribution 

Based  on  Mass 
Median  Drop 

N-2 

40.0 

0.770 

0.805 

N-3 

26.5 

0.690 

0.680 

A-l 

38.0 

0.773 

0.892 

A-2 

30.0 

0.725 

0.732 

A-3 

22.0 

0.632 

0.630 

Essentially  the  same  procedure  was  followed  for  computing  the  overall  capture 
efficiency  of  the  5:1  ellipsoid  of  revolution.  The  method  for  computing  ellipsoid  capture 
efficiency  is  outlined  by  Dorsch,  et  al.  (Ref.  21).  For  this  study,  an  ellipsoid  with  a 
midsection  diameter  of  2  ft  was  chosen;  hence  the  major  diameter  (axis)  was  ten  feet. 
For  the  same  free-stream  condition,  the  computed  capture  efficiencies  of  the  ellipsoid 
are  shown  in  the  following  table: 
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5:1  ELLIPSOID  (10-FT  MAJOR  AXIS) 

Capture  Efficiency,  Overall,  Eco 


Distribution 

microns 

Based  on  Distribution 

Based  on  Mass 
Median  Drop 

N-2 

26.5 

0.107 

0.079 

N-3 

40.0 

0.180 

0.160  ' 

A-l 

38.0 

0.156 

0.140 

A- 2 

30.0 

0.1 19 

0.110 

A- 3 

22.0 

0.072 

0.055 

From  the  data 

in  the  previously 

shown  tabulations,  the  following 

conclusions  can 

be  drawn.  First,  with  only  one  exception,  the  calculated  overall  capture  efficiencies  for 
both  cylinders  and  the  ellipsoid  of  revolution  which  were  based  on  droplet-size  distributions 
are  approximately  equal  to  capture  efficiencies  based  on  the  mass  median  diameter  of 
each  distribution.  The  exception  is  for  the  distribution  N-3  for  both  the  1/2-in.  and  3-in. 
circular  cylinders.  Second,  in  general,  as  the  mass  median  diameter  is  increased,  the  capture 
efficiency  increases,  and  as  the  cylinder  diameter  increases,  the  capture  efficiency  decreases. 
Third,  as  stated,  for  the  ellipsoid  of  revolution  the  capture  efficiency  based  on  the 
droplet-size  distribution  compares  reasonably  close  to  the  capture  efficiency  based  on  the 
mass  median  drop  diameter.  This  agrees  with  the  Langmuir  predictions  (Ref.  26).  It  is 
pertinent  to  point  out  that  Lewis  and  Ruggeri  (Ref.  27)  provide  experimental  data  on 
capture  efficiency  for  four  different  bodies  of  revolution  which  they  compared  with  the 
theoretical  prediction.  They  point  out  that,  in  general,  fairly  good  agreement  exists  between 
experimental  and  theoretical  impingement  characteristics.  Fourth,  the  most  important  point 
for  this  study  is  that  the  results  of  the  capture  efficiency  calculation  showed  that  the 
capture  efficiency  of  the  bodies  was  about  the  same  for  drop  distributions  produced  in 
the  icing  facility  as  for  natural  cloud  icing  conditions  when  the  mass  median  diameters 
for  these  distributions  were  approximately  equal.  However,  the  capture  efficiencies  based 
on  the  test  cell  distributions  were  slightly  greater  than  those  based  on  natural  distributions 
but  not  significantly  so.  Thus,  the  capture  efficiencies  of  the  engine  inlet  components 
under  simulated  icing  conditions  should  be  approximately  equal  to  capture  efficiencies 
which  occur  during  flight.  This  means  that  if  the  ground  test  icing  facility  provides  a 
spray  cloud  with  the  correct  mass  median  drop  size  and  of  the  correct  liquid  water  content, 
then  the  absolute  rates  of  impingement  of  subcooled  water  on  the  engine  inlet  components 
should  correspond  very  nearly  to  in-flight  water  impingement  rates.  Thus,  it  is  possible 
to  provide  very  good  icing  simulation,  in  principle,  in  the  ground  test  facility  for  turbofan 
and  turbojet  engines  of  the  type  which  are  typically  pod  mounted. 

An  important  point  to  resolve  is  the  degree  of  approximation  which  is  acceptable 
when  providing  simulation  of  the  icing  conditions  in  the  ground  test  facility.  That  is, 
how  accurately  in  practice  must  the  free-stream  liquid  water  content,  humidity,  and  other 
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parameters  be  simulated  in  the  icing  facility  in  order  to  achieve  good  simulation,  and 
hence,  accurate  experimental  evaluation  of  anti-icing  power  requirements.  In  the  following 
sections  of  this  report,  the  sensitivity  of  anti-icing  power  requirement  of  a  surface  to 
variations  in  stream  liquid  water  content  and  humidity  is  analyzed.  This  analysis  is  then 
used  to  define  the  kind  of  accuracy  required  in  ground  testing  for  providing  these 
parameters  in  order  to  obtain  good  icing  simulation. 

6.4  HEATING  RATE  REQUIRED  TO  MAINTAIN  AN  ANTI-ICING 

CONDITION  ON  A  SURFACE  IN  AN  ICING  ENVIRONMENT 

The  purpose  of  this  section  is  to  analyze  the  effect  of  humidity ,  liquid  water  content, 
and  liquid  capture  rate  on  the  heating  rate  required  to  maintain  an  anti-icing  condition 
on  engine  components.  The  analysis  parallels  the  analysis  given  in  Ref.  28.  The  degree 
of  importance  of  duplicating  in  the  test  cell  the  values  assumed  by  these  variables  under 
actual  flight  conditions  is  discussed.  The  flow  over  the  component  surface  is  assumed 
to  have  been  accelerated  or  diffused  from  some  prescribed  free-stream  conditions.  Thus, 
the  fact  that  the  Mach  number  of  the  flow  entering  an  engine  compressor  is  generally 
different  from  the  flight  Mach  number  can  be  recognized.  Flow  from  the  free  stream 
to  a  component  surface  is  assumed  to  occur  isentropically  and  at  constant  specific  humidity. 
The  validity  of  the  latter  assumption  can  be  inferred  from  calculations  made  with  the 
computer  program  discussed  in  Section  IV.  For  prescribed  free-stream  conditions,  the 
conditions  at  the  edge  of  the  gas  boundary  layer  on  the  surface  are  determined  by  the 
local  static  pressure  or  Mach  number.  Since  some  of  the  droplets  cannot  follow  streamlines, 
a  certain  portion  of  them  will  impinge  on  the  surface  and  form  a  liquid  film.  The  rate 
of  impingement  depends  on  the  liquid  water  content,  the  surface  capture  efficiency,  and 
the  flow  conditions.  It  is  assumed  that  the  surface  is  maintained  at  the  freezing  point 
(492°R)  and  that  enough  heat  is  supplied  to  prevent  any  ice  formation.  Vaporization  from 
the  liquid  film  into  the  adjacent  gas  boundary  layer  will  occur. 

In  the  absence  of  radiation  heat  transfer,  a  heat  balance  at  the  surface  yeilds 
a"  =  h(T  -  T  )  +  w  (h,  \  +  w  C  h  -  T  ^  -  w  V2  /2 

*w  '  w  aw'  vlil  c  cl  w  c  )  coo 

v  8/  k  oo/  (37) 

I  II  III  IV 

where  q'^.  is  the  anti-icing  heating  rate,  h  is  the  gas  side  heat  transfer  coefficient,  wv 
is  the  rate  of  vaporization  from  the  film,  wc  is  the  flux  of  subcooled  liquid  entering 
the  film,  Tw  is  the  wall  or  film  temperature,  Tc^  is  the  subcooled  liquid  temperature, 
and  TavV  is  the  adiabatic  wall  temperature.  In  this  heat  balance,  it  is  tacitly  assumed  that 
the  liquid  in  the  film  is  at  a  uniform  temperature.  A  few  remarks  about  the  individual 
terms  in  Eq.  (31)  are  in  order. 

Term  (l)  is  the  rate  of  heat  loss  from  the  liquid  film  to  the  air  by  convection.  The 
adiabatic  wall  temperature  is  given  by 
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where  r  is  the  recovery  factor  and  Me  is  the  Mach  number  at  the  edge  of  the  gas  boundary 
layer. 

Term  (11)  is  the  rate  of  heat  loss  due  to  evaporation  of  the  liquid  film.  In  terms 
of  the  mass  transfer  coefficient  (kx),  the  rate  of  evaporation  is  approximately 

wv  =  ^x(Pvw-PveyPe  (39> 

where  pVw  is  the  saturation  pressure  of  water  at  the  freezing  point.  By  the  analogy  between 
heat  and  mass  transfer  on  a  flat  plate, 

_h_/Pr\0-67 

*  C  M  \Sc/  (40) 

p 

Term  (III)  is  identified  with  the  sensible  heating  of  the  supercooled  liquid  droplets 
to  the  film  temperature  from  their  temperature  before  impingement.  The  droplets  are 
assumed  to  pass  through  the  gas  boundary  layer  with  no  temperature  or  mass  change. 

Term  (IV)  reflects  the  flux  of  kinetic  energy  contained  in  the  droplets  entering  the 
film.  It  is  assumed  that  the  droplets  are  in  velocity  equilibrium  with  the  gas  at  the  gas 
boundary  layer  edge  and  that  the  droplets  do  not  slow  down  as  they  pass  through  the 
gas  boundary  layer. 


It  is  convenient  to  introduce  the  parameter  (b  =  wcCc/h)  and  rewrite  Eq.  (31)  in 
the  form 


It  can  be  noted  that  the  parameter  on  the  left-hand  side  (Q)  is  the  ratio  of  heating  rate 
required  in  an  icing  environment  to  the  heating  rate  required  to  maintain  the  surface 
temperature  at  Tw  in  a  dry  environment.  The  parameter  (b)  is  approximately  the  ratio 
of  the  sensible  heat  transfer  to  the  impinging  liquid  to  the  heat  loss  from  the  film  by 
convection.  It  reflects  the  effect  of  liquid  water  content  and  the  liquid  capture  rate  on 
the  heating  rate.  In  particular,  b  is  directly  proportional  to  the  liquid  water  content  and| 
surface  capture  efficiency.  The  humidity  is  reflected  in  the  partial  pressure  (pVe). 

Based  on  calculations  made  with  the  computer  program,  it  can  be  inferred  that  the 
supercooled  droplets  are  in  quasi-equilibrium  with  the  gas  phase  locally  along  the  pathl 
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from  the  free  stream  to  the  compressor  inlet  and  icing  surface.  The  process  occurs  at 
substantially  constant  specific  humidity,  and  depending  on  the  engine  inlet  ducting,  the 
flow  may  be  isentropic  to  a  good  approximation. 

Thus,  for  given  free-stream  pressure  (p^)  (or  altitude),  temperature  (T^),  and  velocity 
(VJ  or  Mach  number  (MJ,  the  heating  rate  parameter  (Q)  can  be  determined  as  a  function 
of  local  icing  surface  pressure  (pe)  or  Mach  number  (Me),  the  parameter  (b),  and  the 
humidity.  The  subscript  (e)  denotes  conditions  at  the  outer  edge  of  the  gas  boundary 
layer.  From  the  constant  humidity  assumption, 

Pv  7PC  =  Pv  /Poc 

e  oo 


Figure  26  shows  the  effects  of  boundary  layer  edge  Mach  number  (Me)  and  the 
parameter  (b)  on  the  heating  rate  parameter  (Q)  for  the  high  altitude  test  conditions  of 
the  military  specification  given  in  Table  I.  Calculations  have  shown  that,  for  icing  surfaces 
in  a  turbine  engine,  b  generally  will  not  exceed  2.0.  It  can  be  observed  that  Q  is  relatively 
independent  of  the  free-stream  Mach  number.  If  M„  is  increased  from  0  to  2,  Q  only 
increases  by  a  factor  of  two.  This  is  an  important  result  since  b  is  directly  proportional 
to  the  liquid  water  content  and  has  some  dependence  on  the  droplet  size  distribution. 
If  b  is  approximately  0.5,  a  50-percent  variation  in  b  results  in  about  a  10-percent  variation 
in  Q.  It  can  be  concluded  that  exact  duplication  of  the  liquid  water  content  or  droplet 
size  distribution  is  not  too  important  in  conducting  anti-icing  tests.  It  should  be  pointed 
out  that  this  conclusion  may  not  be  applicable  for  conducting  icing  tests  without  anti-icing 
heat  or  for  testing  de-icing  equipment.  The  rate  of  ice  accretion  on  an  unheated  surface 
is  highly  dependent  on  liquid  water  content,  liquid  capture  rate,  and  the  frequency  of 
operation;  therefore,  the  heating  load  of  a  de-icing  system  is  proportional  to  the  rate 
of  ice  accretion. 

In  Fig.  27,  the  effect  of  humidity  on  the  heating  rate  parameter  (Q)  is  shown.  It 
can  be  noted  that,  for  the  particular  flow  conditions  for  which  the  calculations  were  made, 
the  heating  rate  parameter  is  relatively  insensitive  to  the  humidity.  This  effect  can  be 
attributed  to  the  low  partial  pressure  of  water  vapor  at  the  edge  of  the  boundary  layer 
compared  with  the  partial  pressure  at  the  surface  of  the  liquid  film.  The  partial  pressure 
difference  across  the  boundary  layer,  of  course,  determines  the  rate  of  evaporative  cooling. 
Large  variations  in  the  small  partial  pressure  at  the  boundary  layer  edge  do  not  significantly 
affect  the  partial  pressure  difference. 

It  is  important  to  have  approximately  the  correct  humidity,  however.  Icing  tests  have 
been  conducted  in  which  the  air  was  essentially  saturated  at  the  total  conditions  of  the 
stream  corresponding  to  a  large  degree  of  supersaturation  in  the  free  stream,  which  implies 
that  the  airflow  into  the  engine  compressor  was  supersaturated  in  the  test.  This  procedure 
cannot  produce  the  desired  degree  of  simulation,  especially  for  flight  Mach  numbers  above 
0.5. 


The  effect  of  humidity  on  the  heating  rate  parameter  becomes  more  pronounced 
as  the  air  temperature  is  increased,  although  the  absolute  heating  load  diminishes.  This 
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effect  deserves  detailed  consideration  if  anti-icing  tests  are  to  be  conducted  at  temperatures 
much  above  the  -4°F  specified  in  Table  1. 

<i 

6.5  EVAPORATION  FROM  THE  SPRAYBANK 

In  order  to  prevent  freezing  in  the  water  spray  system,  the  water  must  be  sprayed 
into  the  cold  air  stream  at  a  temperature  exceeding  the  freezing  point.  Thus,  evaporation 
of  the  spray  will  occur.  Net  condensation  may  occur  if  the  air  is  greatly  supersaturated. 
The  extent  of  evaporation  is  governed  principally  by  the  air  temperature  and  humidity 
and  the  temperature  of  the  water  at  injection.  The  initial  rate  of  evaporation  from  a 
droplet  may  be  very  high;  however,  evaporative  cooling  and  convective  heat  transfer  from 
the  droplet  bring  the  droplet  rapidly  to  the  local  wet  bulb  temperature  and  a 
correspondingly  low  evaporation  rate. 

There  are  two  primary  effects  of  spray  evaporation  to  be  considered.  The  first  effect 
is  humidification  of  the  air  stream.  As  discussed  in  Section  5.4.6,  it  may  be  desirable 
to  create  the  desired  humidity  by  mixing  steam  with  the  air  upstream  of  the  spraybank. 
Thus  it  is  important ’to  know  how  the  spray  evaporation  affects  humidity.  Based  on 
calculations  made  with,  the  computer  program  discussed  in  Section  V,  it  appears  that  spray 
evaporation  is  not  too  important.  If  the  air  stream  is  nearly  dry,  evaporation  from  the 
spray  may  more  than  double  the  humidity.  However,  since  the  evaporative  cooling 
contribution  to  the  heating  load  for  anti-icing  is  not  highly  dependent  on  humidity  for 
nearly  dry  air,  a  relatively  large  change  in  humidity  will  not  seriously  affect  the  anti-icing 
heating  rate.  On  the  other  hand,  if  the  humidity  is  large,  the  percent  change  in  humidity 
resulting  from  the  spray  evaporation  is  small.  For  high  humidity  condition,  the  evaporative 
cooling  contribution  to  the  anti-icing  heat  load  is  small;  thus,  the  additional  humidification 
from  the  spray  should  have  but  a  small  influence  on  the  anti-icing  heating  rate.  The  second 
effect  is  the  change  in  the  liquid  water  content  due  to  spray  evaporation.  The  average 
liquid  water  content  in  an  icing  tunnel  is  usually  determined  from  measurements  of  the 
air  and  water  flow  rates.  Results  of  calculations  made  with  the  computer  program  show 
that  the  liquid  water  content  may  decrease  by  about  10  percent  because  of  evaporation 
for  liquid  water  contents  on  the  order  of  1  gm/m3.  If  the  spray  vaporization  rate  can 
be  calculated  accurately  for  a  particular  system,  then  the  effect  of  evaporation  can  be 
accounted  for  by  increasing  the  water  flow  from  the  sprays.  In  view  of  the  relative 
insensitivity  of  the  anti-icing  heating  rate  on  liquid  water  content  as  discussed  in  the 
previous  section,  this  effect  is  probably  negligible.  An  upper  bound  on  the  evaporation 
can  be  obtained  by  determining  the  amount  of  liquid  necessary  to  saturate  the  air  stream 
at  the  stagnation  conditions. 


SECTION  VII 
CONCLUDING  REMARKS 

A  study  of  icing  simulation  in  altitude  test  cells  was  made.  A  mathematical  model 
for  the  flow  in  a  typical  icing  test  cell  was  developed,  and  the  governing  equations  were 
programmed  for  computer  solution.  A  parametric  study  was  performed  to  determine  the 
effects  of  test  cell  inlet  and  water  spray  conditions  on  the  thermodynamic  and  kinetic 
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state  of  the  flow  at  the  test  section  or  test  article.  The  importance  of  correctly  simulating 
droplet  size  distribution,  mean  effective  droplet  diameter,  liquid  water  content,  and 
humidity  was  investigated  analytically. 

The  results  of  the  parametric  study  made  with  the  mathematical  model  of  the  icing 
test  cell  indicate  that  simulated  flight  conditions  at  the  face  of  an  engine  compressor 
can  be  obtained  through  proper  selection  of  the  test  cell  inlet  conditions  and  water  spray 
settings.  The  model  can  be  used  to  determine  the  required  test  cell  inlet  conditions  for 
proper  simulation.  The  thermodynamic  state  at  the  engine  was  found  to  be  very  sensitive 
to  inlet  air  temperature,  pressure,  and  humidity  and  insensitive  to  the  inlet  water 
temperatures,  droplet  size  distribution  and  velocity,  and  liquid  water  content.  The  velocity 
difference  between  the  inlet  airflow  and  the  water  from  the  spray  system  was  relatively 
important,  however. 

In  general,  the  total  pressure  and  temperature  at  the  test  cell  inlet  should  be  set 
equal  to  the  total  pressure  and  temperature  of  the  free  flight  condition  being  simulated. 
While  there  is  some  humidification  from  the  water  spray  system,  this  influence  can  usually 
be  ignored.  The  test  cell  inlet  specific  humidity  should  be  set  equal  to  the  specific  humidity 
corresponding  to  saturation  at  the  static  pressure  (altitude)  and  static  temperature  being 
simulated.  The  test  article  should  be  located  sufficiently  far  downstream  of  the  spray 
section  to  ensure  that  the  required  near-equilibrium  thermal  state  is  obtained  in  the  test 
section.  The  required  distance  can  be  calculated  using  the  mathematical  model  and 
computer  program  developed  in  this  study.  For  the  system  considered  in  this  study,  the 
distance  is  approximately  20  ft. 

It  is  very  difficult  to  make  accurate  humidity  measurements  under  conditions 
prevailing  in  the  test  cell  under  icing  conditions.  This  difficulty  can  be  circumvented  by 
mixing  superheated  steam  in  correct  proportion  with  dry  air  entering  the  test  cell.  The 
mixing  process  should  occur  well  upstream  of  the  water  spray  section.  Thus,  the  desired 
test  cell  humidity  can  be  created  and  direct  measurement  with  its  attendant  difficulty 
and  uncertainty  is  avoided. 

An  analysis  of  the  heating  rate  required  to  maintain  an  anti-icing  condition  on  a 
surface  was  made.  From  the  analysis,  it  was  determined  that  the  heating  rate  is  relatively 
insensitive  to  droplet  size  distribution  and  liquid  water  content  and  only  moderately 
sensitive  to  humidity.  It  can  be  concluded  that  exact  duplication  of  these  variables  is 
not  required  for  adequate  testing  of  an  anti-icing  system.  Icing  tests  in  the  absence  of 
anti-icing  heat  and  the  testing  of  de-icing  systems  will  generally  require  better  simulation 
of  the  environment  than  testing  under  anti-icing  conditions. 

The  military  specifications  for  conducting  icing  tests  were  found  to  be  inexplicit. 
There  is  ambiguity  associated  with  the  temperature  and  mean  effective  droplet  size  and 
complete  absence  of  a  humidity  specification.  A  wide  range  of  test  conditions  is  possible 
within  the  framework  of  the  specifications  as  formulated.  In  order  to  be  definitive,  icing 
test  specifications  should  include  the  conditions  in  the  atmosphere  to  be  simulated  including 
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the  ambient  temperature  and  pressure,  the  humidity,  the  liquid  water  content,  and  the 
mass  median  or  mass  mean  droplet  diameter.  I 

♦ 

From  a  simplified  analysis  of  droplet  capture  efficiency  on  bodies  which  represent 
various  turbine  engine  components,  it  was  found  that  the  overall  capture  efficiencies  are 
relatively  insensitive  to  droplet  distribution  when  the  majority  of  drops  are  greater  than 
five  microns  in  diameter.  The  capture  efficiency  (scoop  factor)  for  the  engine  was  found 
to  be  near  unity.  , 

It  appears  that  current  capability  to  provide  icing  environments  for  turbine  engine  , 
testing  in  altitude  test  cells  is  quite  good.  These  facilities  can  provide  the  correct  spray  i 
cloud,  that  is,  one  with  the  proper  mean  effective  drop  diameter  and  a  sufficiently  correct  1 
drop  size  distribution,  at  the  required  thermal. state.  In  principle,  control  over  the  pertinent  1 
flow  variables  such  as  humidity,  temperature,  pressure,  liquid  water  content  and  velocity  ' 
is  excellent. 

This  ability  to  define  and  control  the  simulated  environment  gives  altitude  test  cells 
a  distinct  advantage  over  testing  in  natural  icing  environments  or  artificial  environments 
created  by  tanker  aircraft. 
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Fig.  1  Typical  Propulsion  Engine  Altitude  Test  Cell  Configuration 
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Liquid  Water  Content,  gm/m' 


1.  Pressure  Altitude  Range,  Sea  Level  to  22,000  ft 

2.  Maximum  Vertical  Extent,  6,500  ft 

3.  Horizontal  Extent,  Standard  Distance  of 
17.4  nautical  miles 


Mean  Effective  Drop  Diameter,  microns 
Fig.  2  Liquid  Water  Content  versus  Mean  Effective  Drop  Diameter,  Stratiform  Clouds 
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Liquid  Water  Content,  gm/m 
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1.  Pressure  Altitude  Range,  4,000  to  22,000  ft 

2.  Horizontal  Extent,  Standard  Distance  of 
2.6  nautical  miles 
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Mean  Effective  Drop  Diameter,  microns 

Fig.  4  Liquid  Water  Content  versus  Mean  Effective  Drop  Diameter,  Cumuliform  Clouds 
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Ambient  Temperature, 
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0  8  16  24  30  x  103 

Pressure  Altitude,  ft 

Fig.  5  Static  Temperature  versus  Static  Pressure  for  Icing  Conditions  in  Cumuliform  Clouds 
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a.  Inlet  Mach  Number  Less  Than  Flight  Mach  Number 


b.  Inlet  Mach  Number  Greater  Than  Flight  Mach  Number 
Fig.  6  Schematic  Showing  Possible  Streamtube  Configurations  for  Turbofan  Icing  Conditions 
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Ratio  of  Water  Loading  Factor  in  Test  Facility  to 
Free-Stream  Loading  Factor,  r,  f  /f 
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Fig.  7  Variation  of  Loading  Factor  Ratio  with  Flight  Mach  Number  and 
Engine  Compressor  Face  Mach  Number 
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Fjg.  8  Simplified  Schematic  of  Duct  Geometry  for  Analysis 
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Droplet  Temperature  Prior  to  Freezing 
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Fig.  9  Supercooled  Droplet  Freezing  Temperature  as  a  Function  of 
Droplet  Diameter 


'ty> 

I© 


r 


x/R 

r/R 

pg 

Tg 

Mg 

0 

1. 575 

1355.0 

486.00 

0.074 

0.064 

1.492 

1353.7 

486. 98 

0.083 

0.128 

1.419 

1352.2 

485.85 

0.092 

0.270 

1.288 

1348.4 

485.48 

0.  112 

0.339 

1.238 

1346.3 

485.27 

0.121 

0.524 

1.129 

1340.0 

484.63 

0.  146 

0.593 

1.099 

1337.6 

484.38 

0.  155 

0.708 

1.059 

1333.9 

484.00 

0.  167 

0.916 

1.011 

1328.5 

483 . 44 

0.  184 

1 . 10 

1.000 

1327.2 

483 . 31 

0. 188 

4 

4 

4 

4 

4 

10.00 

1.00 

1327.2 

483.31 

0.  188 

Note  : 

Computer  program  requires  geometry 
input  as  dimensional  quantities, 
that  is,  x  and  r  are  input  in  units 
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Fig.  10  Duct  Geometry  Used  in  the  Parametric  Study  of  This  Report 
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Mass  Fraction  of  Liquid  Water  in  Given  Sample 


Gas  Temperature  (Air  +  Water  Vapor), 
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bD 

H  540 


Calculation  Made  with  Droplet  Distribution  A-3 

Calculation  Made  with  Uniform  Distribution, 
d  «=  25  microns 


Both  distributions  have  same  liquid  water 
content  and  mass  median  droplet  diameter . 


Axial  Distance,  x,  ft 

Fig.  12  Gas  Temperature  versus  Axial  Location:  The  Effect  of  Typical  Drop 
Distribution  Compared  with  Uniform  Drop  Distribution  Having  the 
Same  Liquid  Water  Content  and  Mass  Median  Diameter 


Axial  Distance,  x,  ft 


Fig.  13  Droplet  Temperature  versus  Axial  Location:  The  Effect  of  Droplet 
Size  on  Droplet  Temperature 
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Fig.  14  Droplet  Temperature  versus  Axial  Location:  Single  Size  Drop,  25-micron  Diameter 
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Duct  Inlet  Conditions 


Temperature , 
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Fig.  16  Temperature  versus  Axial  Location:  A  Comparison  between  Water  Droplet 
Temperature  and  Gas  Static  Temperature  for  Two  Different  Initial  Water 
Spray  Temperatures 
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Gas  Temperature 
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Fig.  17  Gas  Static  Temperature  versus  Axial  Location:  The  Effect  of  Liquid 
Water  Content  on  Gas  Static  Temperature 
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Fig.  19  Temperature  Difference  between  Gas  and  Water  Droplets  versus  Axial  Location: 
The  Effect  of  Droplet  Diameter  on  Temperature  Difference 
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Fig.  20  Temperature  Difference  between  Gas  and  Water  Droplets  versus  Axial 
Location:  The  Effect  of  Duct  Inlet  Air  Specific  Humidity 
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Temperature  Difference,  AT 
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Fig.  21  Temperature  Difference  between  Gas  and  Water  Droplets  versus  Axial 
Location:  The  Effect  of  the  Initial  Velocity  Difference  between  Phases 
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a.  Capture  Process  for  Droplets  Not  Influenced  by  Streamlines,  Idealized  Process 


Tangent  to  Body  -  AdT 


b.  Capture  Process  for  Droplets  Influenced  by  Streamlines,  Actual  Process 
Fig.  22  Schematics  Illustrating  Droplet  Capture:  The  Impingement  Process 
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Number  Frequency  of  Occurrence  per  Sample,  N. /N 
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Fig.  24  Droplet  Size  Class  versus  Number  Frequency  of  Occurrence:  Comparison 
of  Artificial  and  Natural  Distributions 
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Number  Frequency  of  Occurrence  per  Sample, 


Fig.  25  Droplet  Size  Class  versus  Number  Frequency  of  Occurrence:  Comparison 
of  Natural  Distribution  to  Langmuir’s  Distributions 


Non-Dimensional  Heating  Rate,  Q  =  q"/h  (T 


Fig.  26  Sensitivity  of  Nondimensional  Heating  Rate  (Q)  to  Sensible  Heat  Parameter  (b) 
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TABLE  I 

ICING  CONDITIONS  SPECIFIED  IN  THE  MILITARY 
SPECIFICATIONS  FOR  TURBOJET  AND  TURBOFAN 

TESTING 


I.  SEA-LEVEL  ANTI-ICING  CONDITIONS 

ATTRIBUTE 
Liquid  Water  Content 
Atmosphere  Air  Temperature 
Flight  Velocity 
Altitude 

Mean  Effective  Drop 
Diameter 

II.  ALTITUDE  ANTI -ICING  CONDITIONS 

ATTRIBUTE  CONDITION  I 

Liquid  Water  Content  0.5  gm/m^ 

Inlet  Air  Temperature  -4°F  (-20°C) 

Flight  Velocity  (Mach  0.32 

Number) 

Altitude  20.  000  ft 

Mean  Effective  Drop  15  microns 

Diameter 


CONDITION  I 
1  gm/m^ 
-4°F  (- 20°C) 
Static 
Sea  Level 
15  microns 


CONDITION  II 

2  gm/m^ 

+  23°F  (-5.  0°C) 
Static 
Sea  Level 
25  microns 


CONDITION  II 
0. 5  g m/m2 
-4°F  (-20°C) 
0.  71 

20.  000  ft 
15  microns 
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APPENDIX  III 
COMPUTER  PROGRAM 

A  listing  of  the  computer  program  used  to  integrate  the  equations  governing  flow 
in  a  typical  icing  test  facility  are  given  in  this  Appendix  along  with  a  sample  calculation 
made  with  the  program.  The  capabilities  of  the  program  have  been  discussed  in  Section 
IV  of  the  text. 

INPUT  AND  OUTPUT 

The  English  system  of  units  is  used  in  the  program  with  lengths,  including  droplet 
diameter,  in  feet,  pressure  in  pounds  of  force  per  square  foot,  temperatures  in  degrees 
Rankine,  mass  in  pounds,  and  time  in  seconds.  Almost  all  input  data  are  punched  on 
cards  and  read  into  the  program.  An  exception  to  this  is  the  initial  integration  step  size 
DX  in  feet,  which  is  specified  in  the  main  program.  The  procedure  for  preparation  of 
input  cards  is  illustrated  below.  The  input  data  required  are  as  follows: 

1.  A  table  of  duct  radius  versus  axial  position 

2.  Gas  properties 

a.  Static  temperature 

b.  Static  pressure 

c.  Mass  fraction  of  water  vapor 

d.  Velocity 

3.  Total  number  of  injection  stations 

4.  Water  properties  for  each  injection  station 

a.  Temperature 

b.  Loading  fraction  (mass  rate  of  water  injection  divided  by 
mass  flow  rate  of  dry  air  through  the  duct) 

c.  Droplet  diameter 

d.  Velocity 

The  input  data  cards  are  prepared  as  follows: 

Card  1 

Cols.  1-2  Number  of  points  describing  duct  geometry  (NXY) 

A  constant  area  duct  can  be  treated  by  specifying  three  points  along  the  duct, 
including  the  beginning  and  end  points. 
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Cards  2  to  NXY  +  2 

These  cards  specify  the  duct  geometry. 

Cols.  1-12  Axial  location 

Cols.  13-24  Corresponding  duct  radius 

Card  NXY  +  3 

Cols.  1-80  Title  identification 

Card  NXY  +  4 

Cols.  1-2  Number  of  injection  stations,  NSTA 

Card  NXY  +  5 

Cols.  1-10  Axial  location  of  first  injection  station 

Cols.  11-20  Axial  location  of  second  injection  station 

Cols.  1  +  10(1-1) 

through  101  Axial  location  of  Ith  injection  station 

A  fictitious  injection  station  at  the  end  of  the  duct  must  be  indicated  on  this  card. 
For  example,  if  only  one  injection  station  is  used,  the  location  of  the  end  of  the  duct 
must  be  entered  in  Cols.  1 1-20. 

Card  NXY  +  6 

Cols.  1-10  Mass  fraction  of  water  vapor 

Cols.  1  1-20  Gas  velocity 

Cols.  21-30  Gas  temperature 

Cols.  31-40  Pressure 

Card  NXY  +  7 

Cols.  1-10  Mass  flow  rate  of  water  from  first  injection  station 

divided  by  mass  flow  rate  of  dry  air  through  duct 

Cols.  11-20  Mass  flow  rate  of  water  from  second  injection  station 

divided  by  mass  flow  rate  of  dry  air  through  duct 
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Cols.  1  +  10(M) 

through  101  Mass  flow  rate  of  water  from  1th  injection  station 
divided  by  mass  flow  rate  of  dry  air  througn  duct 

Card  NXY  +  8 

Cols.  1-10  Water  injection  velocity  for  first  injection  station 

Cols.  11-20  Water  injection  velocity  for  second  injection  station 

Cols.  1  +  10(1-1) 

through  101  Water  injection  velocity  for  Ith  injection  station 

Card  NXY  +  9 

Cols.  1-10  Temperature  of  water  injected  at  first  injection  station 

Cols.  11-20  Temperature  of  water  injected  at  second  injection  station 

Cols.  1  +  10(1-1) 

through  101  Temperature  of  water  injected  at  Ith  injection  station 

Card  NXY  +  10 

Cols.  1-10  Droplet  diameter  of  water  injected  at  first  injection  station 

Cols.  1 1-20  Droplet  diameter  of  water  injected  at  second  injection  station 

Cols.  1  +  10(1-1) 

through  101  Droplet  diameter  of  water  injected  at  Ith  injection  station 

In  cards  NXY  +  7  to  NXY  +  10,  it  is  necessary  to  make  entries  for  only  the  number 
of  injection  stations  being  considered. 

It  has  been  found  that  the  most  common  cause  of  program  failure  is  use  of  an 
improper  integration  step  size.  No  criterion  has  been  developed  for  step  size  selection; 
however,  through  experience  with  the  program  the  user  may  develop  a  "feel"  for  selecting 
a  satisfactory  step  size.  Difficulties  are  to  be  expected  in  getting  solutions  for  initial 
conditions  corresponding  to  an  extremely  high  degree  of  nonequilibrium. 

The  program  output  includes  the  following:  axial  location  (X),  mass  fraction  of  water 
vapor  (CV),  area  ratio  (A/AO),  gas  velocity  (VG),  condensed  phase  velocity  (VL), 
condensed  phase  loading  fraction  (FL),  gas  temperature  (TG),  condensed  phase  temperature 
(TS),  static  pressure  (P),  and  droplet  diameter  (D).  The  first  line  of  values  of  FL,  VL, 
TS,  and  D  at  a  particular  value  of  X  are  values  corresponding  to  water  that  originated 
at  the  first  injection  station.  The  second  line  of  values  of  FL,  VL,  TS,  and  D  corresponds 
to  the  water  from  the  second  injection  station,  and  so  on. 
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PROGRAM  LISTING 


The  program  listing  is  presented  on  the  following  pages. 


C  MAIN  PROGRAM  for  spray  cooler  and  dryer  problem 
C  IF  This  PROGRAM  should  STOP  BEFORE  THE  REQUIRED  LIMIT  ON  X  IS  MET 
C  KITH  A  PRINTOUT  OF  'TEmP-FI-  IS  DOT  OF  RANGE  TG«XXX  ,  TS»XXX  *AT  Tm£ 

C  END  OF  THE  OUTPUT.  -E  SrQULO  INCREASE  THE  AMOUNT  OF  RATER  TO  BE _ 

C  injected  (FL). 

X - 

c 

X _ 

c 

X _ - _ 

c 

x _ 

DIMENSION  mol ( 10) «CPL < 10) 

_ PIMENSI0N_FL(1U)  *VL(10)  *TS(1())  .0(10)  tSTA(ll) _ 

DIMENSION  r(Rfl) 

_ DIMENSION  YP (60) _ 

UlMENSION  huFF(iono) 

_ COMMON  BA _ 

common  HUL. VR.GR.GC*RViCJ*RNCA0*B(*) . 1ST A* I SET *KKK » 

_ 1 R QOnT  t  nST  A « I OUM _ 

COMMON  /  ICFE/CPI iROI i SUMJ { 10)  t DELE (10) »SUM9 (10) » NF ( 1 0 )  , NG ( 1 0 ) 

_ COMMON  /XYn/  XS  <991. »YS  ( V9 1-J.N.X  Y _ 

COmmON/LO*/  MlCM 

- aiUMfla-ZutixA/  - 

COMMON  /HSlJR/  R2 

_ CUMJjQN  ./HOYS/  ALP  ( Lftli  JECT  iiNlT _ _ _ 

CALL  LONITS 

_ LCR«29 _ _ _ _ _ 

LPR-21 


READ  (LCR • 1 02 ) NxY 

_ hXAU tLCRj-104>  IXS1X1  ilSIll  _ _ 

10*  format  (2E12.U) 

_ R2»  YS  ( 1 )  +YS  ( 1 ) _ ; _ . 

NXY»NXY-2 

_ 2Q0  h-M* J _ , - i _ 

BA  ■  0.05 

_ M1CH»U _  •  _ ; _ _ 

mEA0(LCR*51015> ALP 

51015  F ORMAT ( 20 A*) _ 

REAO(LC«tlC2)NSTA 

_ NS1»  NST  A*  1 _ t _ 

READ (LCR* 101 ) < ST A < I ) * I"1 *NS 1 ) 

_ HFADIi  CR.lnf))CV.VG*TG*P _ 

READ (LCR* 1 01 ) (FL ( II » 1*1 »NS1 A) 

_ REAP ILCRt 101)  < VL  < I? >1*1 >nSTAJ _ 

READ (LCR* 101) (TS(I) *I»1.nSTA) 

- READ  (LCR  *101)  (t)  ( I )  *  I  ■  1 .  NST  A) _ 

100  FORMAT (*E10. 2) 

101  FORMAT  IBEIO.O) _ _ _ _ _ 

102  FORMAT*  12) 

_ 103  FORMAT  1 15H1  BEOIN  STATION  *I4.12X*16H  RUN  INFORMATION  *20A4) 

CJ«T78.Q 

_ frO  32.2 _ 

1SET  ■  9 

_ 6d»29.0 _ 

VN*lB.O 

_ RV»  1.986*  T7H.0 _ _ 

DO  20  I»l*10 
NO ( I ) *0 
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NF ( 1 ) »6 

_ £C>i.  ( I )  «feg.»OQ _ 

CPL(I>*1.00 

_ _ SUH9  (  I )  »0 . 0 _ _ _ _ 

2l)  SuMj l I ) *0 .00 

_ "OI»ROL(l)»U.9lT(jQ _ 

cpi»cplu>«o.*«soo 

_ ThEF  ■  S40.0 _ 

1ST  4=1 

_ nEQ»b _ _ _ 

C  The  ARRAYS  USED  IN  DIFFE  AWE  N0»  SET  OP 

Y  ( 1 )  *CV _ : _ 1 _ 

Y  <2)»VG 

_ Y  (3)aTG _ _ _ _ _ 

Y (4) aP 

_ Y (5) aFL  ( 1  I _ 

Y(6)«WL<1) 

miiiiiii _ 

Y (81-  0(1) 

_ aNCAQa  Y(2)«(1.0-Y(111»Y(4)/IY(31»RV(  (l.-YdD/GW  ♦Y<l)/VtO) 

x«o.n 

_ »HlTE(LPK.lu3)  ISTAtALP _ _ _ 

INI TaO 

"3  continue — This  card  establishes  initial 

_ px».oool'*"^  integration  step  size.  _ 

K.OUNT  aO 

_ LL"IT-99999 _ 

00  1  lal.LIMIT 

_ CALL  DIFFE I X , Y , YP, PX ,KKK» NEQi U _ 

IF(X.GE.STA(I5TA*1> )G0  TO  2 

J _ CONTINUE _ _ _ _ _ 

JJJa  JJJ*1 

_2 _ CONTINUE _ ; _ 

C  wE  MILL  NOW  PLOT  THE  STATION  JUST  FINISHED 

_ NP»KOUNT _ i _ 

ISTAa  I$TA*1 

_ I  SET  ■  9 _ 

IF ( 1ST  A  .GT.  NSTA1G0  TO  5 

_ WHITE  (LPH»103)  ISTA.ALP  _ 

C  SET  UP  ARRAYS  FOR  DIFFE 

_ y(neo*i>«fl(ISTai _ _ _ _ _ : _ 

Y  (NEQV)  »VL  ( ISTA) 

_ Y(NEu*3)aTS(ISTA) _ 

Y(N£0*4)aD(ISTA) 

_ N£Q  aNEQ*4  .  _ _ _ 

KOUNT  ■  0 

_ GO  To  3  _ ; _ ; _ 

S  CONTINUE 

_ GO  TO  200  _ _ _ 

222  STOP 

_ END _ 

SUBROUTINE  OIFFE ( X  * Y  ,  YP » OX *KK ♦ NEOi I) 

_ DIMENSION  Y (50) >YP(50>  tZ (501 »2P<50) .ZN(SO) _ 

COMMON  HA 

_ COMMON  ROL  t  VYu  GW  >GC>RVtCJ>  WNC A  0  »  B ( 4 ) 1 1 ST A t ISET  «  KKK  > _ 

1K0UNT  *NST A  t IOUM 

_ COMMON  /XYN/  XS(99) t  YS ( 99 ) »NXY _ 

1  COMMON  /doxx/  dxx 

_ COMMON  /BARK/  B ARH (10) _ 

CALL  LUNITS 

_ LPRa21 _ _ 

LCH»29 

_ CALL  MQNITH(X»Y) _ ; _ 

CALL  YFUNC (XtY*YP» 1 ) 

_ DO  1  I  a  1  .  NEQ  _  __  _ _ _ _ 


73 


AE  DC-T  R-73-1 44 


1  A(I)aY(I)*OX«YP(I) 

Xa  X  »  UX _ 

L)X2a.SO*0*DX 

UQ  S  Ja2t9«?q>j _ r_J _ 

CALL  YFUNC U*Z*ZP» J) 

_ Hill _ 

oo  *0  I«1*neu 

zNin»Y(i)»Dxi>«(YPH)*2P(i)i _ _ 

IF  (UAHS(ZN(I)”Z(I))“1.D“0S*UABS(ZN(I)))4*4*3 
_ 3  Aaj _ 

4  AH«J 

_ _ _ _ _ ; _ 

40  continue 

_ if.  ixi5 ,»  (>  *5 _ : _ 

5  CONTINUE 

_ »KITF (LPH.QRl  (ZNH)  * lal.NEUl _ 

»HITF (LPR*90)KK 

OH  FUHMATlim  _ \ _ 

9R  FORMAT (4E20.10) 

_ hlUfi _ : _ 

6  DO  7  I«1.NEU 

_ 7  mjLMZi-U _ _ 

IF  (J  .GE.  3  .AND.  J  .LE.  5) GO  TO  1212 

IF(.I  .1  T.  31KO  TO  ?Q?C _ 

DX«.b»DX 

_ UXXaOX _ 

GO  TO  1212 

202l)  UX«2.Q»DX _ _ 

IF (OX  .GT.  0.1)  UX-0.1 

_ r.n  Til  \9\9 _ \ _ ’ 

1212  CONTINUE 

_ UAAfQZ; _ 

RETURN 

END _ _ _ 

SUBROUTINE  YFUNC<X,Y,YP.KI 

■  DIMENSION  HOLUfll  tCPUlU) _ 

DIMENSION  Y (SO) .YP(SO) iAb(10l *A6 (10) »A7 (10) tAB{10) • AlO ( 1 0 ) • X VS ( 10 

n.XMl  tlO)  .XMFG(IO)  .BARHIIO)  .ORllO)  .XKI10)  .CDUOl _ 

2  *  XHV (10) 

_ COMMON  HA _ • _ 1 _ 

•  COMMON  R0L»V*.G**GC«RV*CJ»*NCA0«B (4) *ISTA*ISET*KKK* 

_ IMlUNIiMSTiiimm. _ • 

COMMON  /RSOR/  R2' 

_ COMMON  znn 1XL  DA _ : _ : _ 

COMMON  /ROYS/  ALP(lO) • JECTt INIT 

_ COMMON  /ICFF/CPI.RQI.SUHJllQ) t OELE ( 1 0 ) . SUMO ( 1 0 ) . NF ( 1 0 ) .NS ( 1 0 1 

COMMON  /BARK/  8ARH 

_ common  /PI  OTT/  xx  1 1 n?>  .rv  11021.  va  ( 1021 .1011021  .pmo?)  . n do.  102) 

C*VL  (10*11)2)  *TS  (10*  102)  *0  (10*102) 

_ COMMON  /AYE/  I _ 

CALL  LUNlTS 

_ '  LPR»21 _ _ _ 

LCR-29 

_ PIE  a3.141A _ 

SUMlaO. 

_ SUR2«Q. _ 

SUM3*0. 

_ SUH4a  0  « _ 

SUMSaO. 

_ CAIL  BQGCF  1RQG  t  Y  1 1 1 1 Y  12 )  *  X.] _ 

00  1  lal.ISTA 

_ Ja  B.4»(I-1) _ . _ 

A5(l)  a  P1E»R0L(I)*Y( J)«»2  /2,0 

_ CALL  SUBl(Y(l)*Y(3).Y(4).YU-l).XVSm.XV»XK(l)  .  Y  ( J)  *DM(1)  *CO(  I)  *B 

1A«M(I)*Y(J-2)*.Y(2)) 

YP(J)a  XK(I>«PIE«Y(J>«»2.0MXV-XVS(I))«VW/<A5(I)«Y(J-2)M1.0-XVS1I 
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1))) 

_ AMin  ypi.j) _ _ _ 

YP ( J-2) = ( (PIE/M.O) *ROG* <Y  < J) **2  ) *CO( I > *DABS ( Y (2) -Y ( J-2)  ) * ( Y (2) - 

lY(J-2»  »-(  Y(J-2)-Y(g|  )*Y(J-2.)*A5(I)*Ab(I)  )/(DH(I)«Y(J-2)> _ 

A7< I)«YP(J-2) 

_ TP(J-J)«  A5(I)«A6(I)*Y(J-(]/DM(I) _ 

AH(I)=YP(J-3) 

_ j>UMl»SuMl«A*U)»A6(I)»Y(J-3)/0H(I) _ 

call  mli ( y ( j-i >  «xhl ( i > ♦ ahfg< i ) «xhv  <in 

YP(J-i)  «  RAKH(I)«  PIE  •  Y(J)»»2  *(Y13)~Y(J-1) )/  I  DM ( 1 ) * Y ( J-2 ) ) - 
CAHL  ( I )  *A5  (  I  )  *Af>  (  I)  /OM  (  I )  ♦  Ab  ( 1 )  *  Ab  ( I )  *XHV  ( I )  /l)M  ( I ) 

_ A  1  U  11)  s  YP  (,)-!  ) _ ; _ I _ 

SUf'2  ■  SUM2  ♦  Y  ( J-2)  «Ad(  I)  «  Y(J-3)*a7(I) 

_ iiiPJ  ■  SUM3  «Ad(I)«<XHL(l>  ♦  Y  LJ“2)  •  «2  /  (  2  «  Q*6C*C J )  ) _ 

SUM*  ■  SUM*  ♦ Y ( J— J ) •  Y ( J-2) *A7 < I )/ (GC*CJ) 

_ 1)0*5  «  SUMS  ♦  Y(J-3)«Al(MI) _ ! _ 

1  CUfTlNUE 

_ yp(1)»-sumimi.o-y  (l) )  _ 

A9=YP(1) 

_ CAIL  AAA3(X«Y(1)  fY  (2)  <A3<AA0tQAA0  .  RUiiJ _ 

CALL  OtrU  (  *06  »Y (I) «Y (21  *  Y (3) tY (*) . aR • SUM2 *  SUMS ♦ SUM* «SUM5* A0,a3» A  I 

Cl  tAl?«Al3»Al*«M«CPV.CPA,A2iA*f AAOiA) _ 

YP<  3)  ■  (All*Al*-Al2«Al3)/(All«(CPV*Y(l)/( 1 ,0-Y ( 1 )  ) *CP A ) - A 1 2« A 0 «R06 

_ i  •  k  :  _ 

Alb  =  YP ( 3 ) 

_ YP(2)  «  (A13-AO«ROG«R»AlS)/All _ 

Alb»YP<2) 

_ CALL _ HLh P (KUG.Y(l) .Y(2) , Y ( 3 ) , R . A2 » A3 ♦ A* » A9 . Al 5 1 A 1 b .DPiA) _ 

V  P  (  *  )  ■  DP 

_ 1PIK.KT.  11  GO  TO  3 _ _ _ 

IF  ( A  .FQ.  O.OIGU  TO  9999 

_ 1FIX.LT. 0.00151)  OQ  TO_9999 _ 

IP (X  .LT.  nA ) GU  TO  3 

_ IF  (X.LT.1.0) HA  »  BA  «  0.05 _ 

IF  (X.GT.I.U)  Ba  ■  HA  ♦  0.2b 

9999  CONTINUE! _ 

CPG»Y(I)*CPV  ♦ ( 1 .DO-Y ( 1 ) ) *CPA 

20  Format  (in  t3F2u.9) _ ; _ 

CPG»CPG*  7  78.00 

_ 6ah6«CPG/  (CP6-R) _ ' _ 

AMACH-OSORT (GAMG*GC*R*Y (3) ) 

_ AHACH«Y(?)/AHACM _ 1 _ 

P0T»1 .00* (GAMG-1  .00) *0 . 500  * AM ACH* AMACH 

_ PU»Y (*) *POT«* (GAH6/ lGAHG-1  .DU)  ) _ ; _ 

T  U»Y (3) *POT 

30  F0hm«T(1IM1STATI0N  ■  »I*»10X>IfeH  RUN  I NFORM AT  ION » 20A* ) _ 

31  FORMAT (I 5H0 IN I T I AL  AREA  ■»GI6.6) 

_ AHEAI»3.1*15926P0»R2 _ 

IF ( INIT  .EO.  0)  WRITE <LPR*3I  >  AREA  I 

_ INIT.INLL+I _ , 

WRITE (LPR.*)  X • AMACH. POt TO 

_ IF ( JECT  .GE.  &*)  WRITE  <LPRt30) I  ST A I  ALP _ 

IF (JECT.GE.5*)  JECTwO 

_ JECT» JECT*2 _ 

ww»  Y (I ) /  (  I  .0-Y (  1  )  ) 

*  FORMAT(»HOX»  «Glfe.b»3X«*H  MGw«Glb.b»3X.*H  P0»fG16.6> 3* i_*M  T0». 

1  61 6 . b ) 

_ WRITE  (LPW.5)  AAOiWW _ 

IF (JECT  .GE.  5*)  WRITE(LPR»30) ISTA«ALP 

_ IF1JECT.GE.S*)  JECTwp _ 

JECT»JECT*2 

5  F ORM AT ( 8H0  A/AQw  «G16.fe»3X.9H  WV/WNC  ».616.6) _ 

WRITE (LPR  »b ) Y ( I ) «Y(2)*Y(3)»Y(*) 

_ I F ( JECT  .GE.  5*)  WRITE (LPR . 3 0)  ISTA.ALP _ 

IF ( JECT.GE.b*)  JECT  wo 
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b  t  UKM  A  T  (  *  HOC  V=  *  o  1 6  •  b  *  3  A  1 4H  VG= • G 1 b • 6 * 3X * 4H  T  G*  *  6 1  b  .  b  i  3*  •  3H  PB*G16,6 

_ 1 _ I _ 

LpO  10  1«  1  *  1 57  A 

_ .JtHtt*  I  !  -  1  ] _ 

L  E  L  7  =  Y ( J) -Y ( J-l  > 

_ .hTTF  (I  PW.Hl  Y  1.1-3)  .Y  I.I-?)  .Y  1.1-1  )  .  Y  1,1)  iHEl  T _ 

IK  ( J  £  C  7  .  Gt  .  54)  «H17£  (1.PP.30)  157  A*ALP 

_ _ It  urn  .gf  .so  jt ctg o _ 

JtCT= JEC7  *2 

_ ri  FQtMA7 (4H0FI  c.Glo.h.  IX .4H  V I  B«GIfe.fe.3X»4H  T  S»  .  B 1  b  .  b  .  3X  .  3H  I7». _ 

lblb«b»lX»bH  n£LT=.E7.J) 

_ «H1TF(I  PR.U)  YP  1.1-3)  .YP  (.1-/1  .YPM-ll.YPI.il _ 

IF ( J  £  C  T  .Gt.  54)  wh 1 7£ ( LPN • 30 )  IS7A.ALP 

_ if  (Jt  cr.GF.s4)  jtxiaii _ ! _ 

jEC  T * JEC7  *2 

_1 0 _ CODLLUtliE _ 

9  tUt-MAT  (5H(lELF  =  iGlb.6.3A.5H  VLP» * G 1 b . 6 t 3X . 5H  T SPs t G 1 6 . 6 * 3X * 4h  OP»* 

- 1  G  I  b  .  b  )  _  _ 

hHI7f(LPh«7)  YP(1)  ,YP(2)".YP(3)  *YP(4) 

_ IEIJFCT  .fit.  54)  4  P  T  T  t  II  PR ■ J  U  J  ISIA _ 

IF  ( Jt  CT  .GE  .54  )  JEC 7*0 

_ JtCTaJf  CT  ♦  / _ _ 

7  FORmaT  (  5h OC vP*  » G 1  b .  6  *  3A  .  br)  VGP»  *  G1  b.b  *  3X  *bH  TGP« . G1 b .6  *  3/ *4H  f>P»* 

1  G  )  f.  .  b  1 

lFIKOUNT  . GE .10  0)  GO  TO  1000 

_ KOUNT  «  KOUNT  .1 _ 

AX (KOUNT)  3| 

- CM  1KQUN7  1  ■  Y(U _ _ _ _ — - - - - 

VG  <KOUNT)«  Y(?) 

_ T.ii..LnmiNl.l.x  Y.3-1) _ 

P  (KOUNT)  »  Y ( 4 ) 

_ Ull_l0Ui— 1» . 1  jl-LSJLA - 

J=  H*4*(I-1) 

_ El _ (T. KOUNT)  rY(J-3) _ 

VL  (It KOUNT)  BY(J-2) 

_ TS  (T. KOUNT)  «Y(J-1) _ 

0  U.KOUiNT)  «  Y(J) 

...1001 _ LUtiT  1 NUE _ 

1000  CONTINUE 

_3 _ COM  I  NIIF _ 

2  t  OWN  a  T  (  1H  .7tlb.(S) 

_ Kt I Uhttt _ 

END 

_ SUHHUUTTNF  MlmllHU.Y] _ 

DIMENSION  YP ( 50 ) 

_ dimension  xus(iii).x<(io).uw()o).ro(ioi _ 

0 1  MENS  I  ON  Y (50) .b(10) iCPL(IU) 

_ common  Kfl.wni  (  l  ii )  .  v  «  .  i;»  .  gc  .  nv  .  r.i.  wnc  an  .  7  (4  )  .  is  Ta 

COMMON  /let  E/tPI  3  Hoi  f  SUMJ  <  10)  *l)Elt  (  10)  *SUM9<  10)  *NF  (10)  tNbUU) 

_ common  /nnxx/  da _ 

COMMON  /AYE/  I 

_ COMMON  /bAKK/  OAKrl(lO) _ 

CALL  lunits 

_ 1  PPePl _ 

LCR*?9 

_ t  nTI  0»1  44  .  Ii _ 

CALL  YFUNC(AtY.YP,2) 

_ UO  10  1  =  1.ISTA _ 

b ( 1 ) shAKHI 1 ) 

_ 10  LUHTtNlJF _ 

00  1000  1=1.1STA 

_ J=b*4*  I  T-  1 1 _ 

A*(Y (J)/. 001093) 

_ lF.tYlJ)  .I  F.  .(JU32BOn  .AND.  YU)  .GT.  .0010 93)  TSF»463.2D0 _ 

1 F ( Y ( J )  .IT. .00109300)7SFB<463. 20*16. 320*OLOG1A> > 

. . <TSF-492.(U _ _ _ _ _ _ _ 
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1MNGII)  -  1)44*55*55 

44  IFIYIJ-1)  .GT .  TSFlGq  Tu  1  OiLQ _ 

55  IFINFIl)  .tO.  I ) 60  TO  1000 

hW  =  l  349.60-. 55  714290_«J  (J-ll _ 

SUMO (1 ) =  SUM9  I I) ♦ Il3.0/t (J) ) *YP( J)*OA*hVV) 

bUM  J 1 1  )  =  SUM J (  I  )  ♦  6  .  L)Q *fi  ( 1  )  ♦  1Y  (  J  - 1  )  -  Y  (3)  )*DX/(HOL(I)*Yi?'*Y(.j)i _ 

•  illTF  (LRK,  IOC)  SUM9  |  I  )  *  SUM  J  1 1)  •  1 

»HlTf  ILPh,  IQO)  YP  ( J )  .MVV.ISIA _ _ 

»«I  T  t:  (LRH,  101  )  Mil) 

lOl  F  QH4AT  (E2u. 1  U  ) _ _  _ 

ioo  format  i2t?o.io*  no) 

_ SuHJ(1I=SUm.m1)-SUM9(1) _ 

»nl TF (LPR.99)  bUMJ(I) 

_ IF  (UatjS  (SQmJI 1 )  ) -  E  N  T  L  ) 4*5* 5 _ 

4  Y  IJ-1  )  *49?. '10 

_ ■-(4(11*) _ _ 

uu  Til  IUOo 

5  NF<1)*I _ 

*OL I  1 ) *HO  I 

loop  CONTINUE _ _ 

99  t  OPMaT  |F  21).  10) 

_ _ IF  LA.  «  oT  .  10U.19)«WlTL(LPK*99)SUHjm _ 

MtT  URN 

_ 1-  Nil _ 

bUl.HUUT I NF  SUMlI  CV  .  TG.P . TS* AVS* AY* AK.O.Om.CD.BH. VL * VG) 

. _ L'l^FNblON  TT  (21 _ 

U  I  MENS  ION  HOLIlO)  ,CPl_!lO) 

_ COMMON  HA _ _ _ . _ 

COMMON  HUL* VH*bN*GC*RV*CJ*«NC«0*H 14) • 1 5 T A  *  I SE T * KMv  • 

_ lAUUN.LlflS.I.AilLlU£! _ 

COFMUN  /ICE  I ,/CPl.HOl.SUHJIIO)  .06l£  I  10)  *SUH9lI0)  *NF ( I  0)  ,NGI 10) 

_  COMMON  /AYE/  1 _ 

CALL  lunits 

_  LPK»?1 _ 

LCH»29 

_ U£j _ 

IF ( T  S-49? , 015*5*4 

4  X  V*  (CV/ V*»)  /  (CV/Xb  ♦ ( I ,-CV ) /Gx)  _ 

A2*-9.U6*  I  t-.669h*TS*  .0639E-4  *TS**2  ♦ • 0927E-7  * T5**3» 1 35? . 3) /TS- 

_ Cl. 4525744)  _ 

*1=  672.0/Tb 

FVTS  *  2117.0  «Al*«5.19«OEXP(Ag) _ 

00  To  6 

5  IF(NFU) zlltLllxI _ _ _ 

6  A=b. 426651400 

_ TMJ)»T5/1.00  _  ___ _ _ _ 

Frt=-2005.IUO 

_ C= 1 . 3B69L1-4 _ 

DO* I .19650-11 

_ uK«2937o0.U0 _ _ _ 

t  *-.004400 

_ L*- .005714800 _ 

M*I 0.**  I (00* ITT  I J) **2.— DM  **2.0) ) -1 ,0 

_ h»C* (TT IJ) **2-0K) /TT  IJ) _ 

U* 10. 0**1 IK*IJ74.11-TT(J) »  ** 1 .25 ) ) 

_ hq*m*h _ _ 

b  =  A ♦ RR/ T  T ( J) *HK 

_ PvTS»21l7.o*OEXP(2.3* IS*E*U)  ) _ 

GO  TO  « 

,7  _ 

TT (J)*TS/l.dO 

_ BR»6. 231200 _ _ 

C*-. 0167700600 

_ DO*. 0  0  0  0120514QO _ 

£*-6.75716900 

_ _ F»((A/TT(J)  )*^434l»Otf»OL06imjm-iC»TtU).«DD*IJ<0>-«»a*E  . 
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f=2.3*F 

_ k>TS8^.7Hih*[iKAH  t  A  I _ 

K  (.OAT  INOF 

_ XVb  -  PVTS/P _ _ 

A  V  F  =(AVS*AV)/2.0 

_ A.-1SA  lit  S±JL _ «  <  1  AtfF)  »fi« _ 

IF  I  =<lS*TG)/2.0 

_ t  PV*  mi.FiF-7  «Ob(,pr  (TFI))/(l.*lb3H./TF!) _ 

F  *NC  =  7.bE-7  *IjSWHT(TF  I)/  1 1  .  *21  f>  •  / 1 F  I  ) 

_ LK _ =  XVF  *E±LK  *J  1  .-XVF  )  •  F  M  N  C _ 

n»H  *  (  .S37b/P>  «  (  TFI/4SI1  .0)  **2.334 

_ mlF  =  (-•  Art/  IHVMF  I  1 _ 

bCF  =F>'/  (H(»F*|'ah) 

_ ifiva  .fq.  vi  )  oin^t  .oonouoi _ 

Nt  *  1 1  AhS  (  Vj- VL  )  *hOF'*‘)/F  A 

_ xH.-.=  2.0*q.h«iiS'jkT  (Hf  )  •  S(.F~ •  •  U.333J _ 

A  A  =  AhN*FaV  (i>*SCF*AM) 

_ I  ■"  a  n«».i _ F(.U1(1  *HUl  1  T  )  /FwO _ 

UJ  =  24.0/Ht‘«<  1.0*.li*Mt**U.tiB/) 

_ IF(  T  F  1  .Gl. 17(10.0  .AND.  T  F  1  .1  T.diiHil.ll)  GO  TO  221 

IF  (  T F I  .Ll.  400.  .ok.  TF  I  .GT.  4bU(J.  )  00  TO  222 

_ UiY£  =. 4.1()(.  4.  lt.7nF-i.ATFT  (O.A7>OF»7*TFT»*?.(i _ 

CPNCF  *  11.231k  ♦  0.1040F-4  •  I F  x  ♦  0.716AE-8  4TFI**2 

_ nil  In  1  <1  S 

221  CPvF  =  ,3TH  «n.l4jMF-j  • T F I  -0.1312t«?  *TFI**2.0 

_ LPr.CF'*  0 .  22  1  4  ♦  U.JS21F-4«TFT-  Q.J77AF-H  4TFI«42 _ 

o()  T'l  102 

_ 2.2.2  »wll  ^ll.Ek.illiXl _ 

101  FoKmaT(2SH  TtMK-FI-lS  OUT  OF  FtANGE  ) 

- Art T  TF  11  PK a  1  2-2.1-I-lii.I.b _ _  _ 

122  F  OHM*  T  (Sfki  [(,=  .t2u.6«bX«4M  Tb=»E20.o  ) 

_ Al OH _ 

102  CONTINUE 

_ FKV  «  O.X7«h.V _ 

fNnC  *  0 . 2b  7  •(0.11b*  b  .  1  /  •CF’nCF  )  *  F'M.VC 

_ m  *  xvF  «FKV  « I ( .- X  vF  )  »FK.vC _ 

CPF«  XVF* (V./AP) *CP VF  ♦ (1 .-XVF) *CPNCF  *G*/XM 

_ PVF  «  CPF  *F»./FK _ _ 

F NU  =  2.  *0.<'*l'b0HT  (HE)  ♦  PVF  **0.3333 

_ ert  «FNU  «F>,  «CPF/  (U*PVF  ) _ _ _ 

Kt  T  UHN 

_ LhdJ _ _ _ _ _ 

SUKKOuT  INF  HOOEE  (HOiitCV*  VGi  A) 

_ (iTmfmston  kil  (10)  . r.Pi  LLuJ _ 

ul  Mt'Nb Ion  a  ( 4 » b  ) 

_ CUhwch  h  ft _ 

COrtViM  HUL  *  v»  « Op  » l;C »  W V  «  C  J *  *NC  AO  « B  ( 4  )  »JSTA 

_ (■'■)>  MI1V  /XYN/  XS  (SIN)  .YS  (HQ)  .NAY _ 

/Abort/  K 2 

_ CAU  LUJV1T.S _ 

LPP=2l 

_ (  CH  =  2<rt _ , _ 

IF  < A.LT.ASI  1 ) )  GO  TO  lo 

_ if  i».i  t.iupii  r,u  Tn  2d _ 

OU  b  I  *  1 ♦ NX V 

_ IF  1 A.GF .AS  1  T ♦ 1  )  . ANP.X.I  E .xs ( 1*2)  )  SO  TO  30 _ 

IF  (X.OT.XS(NXV) )  GO  TO  10 

_ b  CONTINUE _ 

w«ITF (LPK.lOO) X 

_ 100  FOHMATI17H  X  DOT  OF  H ANGF . » . H fe . 7  ) _ 

STOP 

_ 10  (.UNTInUF _ 

1=nXY-1 

_ 00  TO  3(i _ 

20  1  =  1 
30  CONTINUE 
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n  =  I-1 

: _ uo  ji  i* g 1 1 4 _ 

aIa  =1*K-1 

XN=K«N _ 

MKdlil. 

_ k(K.SJ  sYS(KIK) _ 

UU  3]  L*2  « •« 

_ 1FCXS(kmi  .fo.  tj.)  A  (K«l  )  =  L).LiO _ 

IF (AS(FN)  .EU.  0.)  HO  TO  .11 

_ AlK»L)sAS(KN)««(L-H _ 

31  CONTINUE 

_ LtUiiA _ 

1CL»6 

CALL  ChOLES  IA« IPUtlCLtO) _ 

lFdwU  .FU.  -1)  »«ITE  <LMM«  M)0> 

TOO  FUFMxTtSbU  IhF  -A  T^ I A  USD  To  GET  FIT  OF  DUCT  UFQm£THY  IS  SINGULAR 
1) 

_ UU  t)0  1*1.4 _ 

60  b(I)*MI«S) 

a  a  ti  —  h  ( i )  ♦ki2)«a  h  ui  _ 

»aosaau*aau/«2 

_ hUO  S  *NCA0/t(l.-CVI«V0«A4U) _ 

Kt  TuHN 

_ LM! _ 

bUrtHOUT 1NE  ML  I  (Tbt AHI_* AMF  0« AbV) 

_ XML  =TS  -b4<).0 _  _ 

XriFO  ■  -,SbVt>*TS  ♦  .0rt3Vt-4  *TS**2  ♦  U.0927t’-7  *T  S**  s .  0*  1  3*>2. 3 

A  M Vs  AML  ♦  X  HF  fa _ 

KETUKW 

_ LiiO _ 

bUbHOuT  INE  AAAJ  (X*CV,Vb*A3*AAO*OAAO«NO(») 

_ UIwEnSIQN  KULllOl  iCPLtlol _ 

CONHUN  bA 

_ LUMWON  PUL  i  Vl>  i  Q«  iGCiKV  if  Ji  <1 CA  lliH  (4  )  »  IS  T  A _ 

Common  /*Sup/  n>t 

_ AAU  S  b(l)  «  h(2)«A  *b(3)»A««2.U  «  b  I  a ) »X**3 . U _ 

UAAO  S  H<2)  *2.U*  h(j)*A  ♦ 3.0*8 (4) *X**2.0 

UAAOs2.1)U*AttO*u'AAU/><2 _ 

AA0*AA0*A A0/«2 

A3  =  Vb*(l.~CV)*UAaO/(vl.*ll.-CV)*AAu)*,g« _ 

hCTum'i 

_ LliU _ 

hUbPOUT INF  CHOLFS ( A*N»M»MA TSYM) 

_ UIHElNSION  MN«n) _ 

IF  <  A  { 1  *  1 )  .nF  .  l).o)  (.0  To  <*7 

UU  37  J»?.N  _ _ 

IF  1  A ( 0  *  1 > .FU.U.U)  00  TO  37 

_ IF  L1H.=J _ 

00  To  i7 

37  CO N T  TNUE _ 

00  TO  b«J?l 

2  7  uO  &7  mli" _ _ 

IEhPsA  < IFLIPtK) 

_ A  t  IFLIF*f*  )  g.A  I  1j  K ) _ 

AlltKI STEMP 

b7  CONTINUE _ 

47  00  2  J»2*m 

Ml.JItftU.JI/Mhll _ _ _ 

2  CONTINUE 

_ LilL  fe_I  =  2iN _ 

UO  7  J  =  2  * M 

IF (MATSYM.KU.O)OO  TO  49 _ 

IF  (  I-J) 4 V  »  68  »  6  7 

49  IFIJ.GT. 1 ) m  TO  69 _ 

6B  AsJ-l 

_ SUMs^.  O _  _ _  .  _ _ 
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L  u  3  1«=1 

I.IWlOulIJ.  It 


3  VimT  I  »UF. 

<1  I  1  .  .  I  )  =  A  I  I  .  -I  I  - 


uo  rn  7 

1-1 


bu  (1  =  SUM*fl 
i  i  ifu  T  r 


\¥  (ACI.I) .fj.y.n)  bO  TO  3*321 

/  I  _  I  I  —  I  I  I  .  I  ’  —  l-  l  i  ..1  1  /  I  I  I  .  1  I 


oo  TO  7 


7  ConT  (MiF 


ou  b<^  **2.n 


buMSI'  .  II 
I  I  S  I  A  I 


oO  bl  InsLL*'' 

i  1  .  1  i  (1  .  M 1 


51  (.onT  1  i'uF 

II.I.UAll.M. 


52  Cni'llNoF 

i.(i  Tn  l 


54321  Iv  =  -1 


1  A  i  2  t  4  1  3  *  A  1  4  AflOtA) 

•  tia.ivi  iim  l«  i  ■  i  i  l  t  .i  ui  t  l  . .  i 


4  13*  AU«"NOO»  I  0*0*  A  3- A '1*  (  VG*«2 ,  »KOG* To*  A2* AO  ♦  TO* H*  A*  A  0  )  -VG*  (  1  ,-CV 
1**2.*' 


call  Mf  pv  ( ro* *f.v  »c* v*cha  ) 


Hr!  1  U«U 

in  ii  I  i  r 


LRR  s2 1 


Hit  1  rt  (LFR.Jl  )  V6 

~  ~  i  ii  i 


OUkHOOT  lot  nt-.fc  V  (  TG*XrlV  .  CFV  tCP  A) 

liImFnSTON  ►ill!  [lull  CFL  ( 1  h  ) 


COF'mO'4  HA 


IF (  TO  .GE.17O0.  .AND.  Mi  .Lt.  45UU.1G0  TO  222 
IF(  TO  -IT.  Atm-  .i)«.  1«  .IO.  ARUO.loO 


XMV  s  ,43<>4*TG  *.U«  39F-«.*TO**2.0*0.  J92 7E-7* TG** 3 . 0-236. 3 1  6 1 ♦ 1 04 2 . V 
\i  =  .tint  *  .  1 1,7nf  .4*1  I,  «  .  ?7H  1  P-7*  Tli** 


CHA  =  . 231 H*  .lO*Ufc-4*TG*.71f>6t--M  *TG**2.0 
00  To  H>0 
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22?  amv  =  ,3319  «l(i  ♦  .07i9F-JMi»*«2.o-,U4j/3t-7*TG**J-le5.381*1042.9 
_ CFV  S  -3319  *  .  lit  3hF-3«1G-  0  .  1 3  1  2E.- /•  U _ _ 


CPft  =  .2214  ♦  0 .  jS2  1  4*  1 6  -u.  17  76t-H*TG*«2.0 

bU  TO  1QU _ 


221  CAIL  LUM1TS 
_ I  Ph»?1 _ 


LL6=?9 

wallL-U-PHi-^L 


V9  TU66flt(24h  UmP-G  IS  0U7  OF  HANGE  ) 
_ fhITf  (LP6. 1 1 TG _ 


1  OHM A  1  (4huTb=.Fl6.6) 

JjlUt _ 


lou 


CtlAjT  l  6U F. 


t.wl) 

bUnH'  )UT  I  »|E 


OEHP  (FUb.CVt  VGiTG.H.t  A2.A3.  A4  .  A9  .  A  1 5 1  A 1 6  .  LIP  .  A  ) 


u  1  «r  '.S I  On  rtiJUlui  fCPLMOl 

LUMw|)'j  HA 


COMMON  6UL  t V*«b» .GC.6V iCJi WNCAO .H (41  . 1ST* 

HE*  nlU^h^Alfe  ♦  (FOG«lb«A2*  TG«N«A  )  «_A.9-_tG«> 
htTUWiV 

UiL! _ 


seof 

*?3. 


0.0 

0.4168  Z_ 


6.1418/ 

6  .  62500 


0.63333 

JU2iluiiL_ 


5 .20m 33 
_4.6  7SOO-- 


1.8686/ 


4  •  6?S  il  O 


2.QU33J- 


.4»3>l-Ui 


2  •  5  u  0  0  o 
_2j.91fe.hX- 


4.2166/ 


3.333.1J 


3.99167 
3.9 loo? 


4.16667 

4.26333_ 


3.90333 
■  3»90flIL- 


4.56333 


3. 90  oil 

3.90inj 


5.4)66/ 

5.63333 


3.9090 

JLfVUHU. 


6 .6666  / 

7.5UH1U) 


3.91)1111 

3 . 9  0 1)  11_ 


8.33  •»  3  J 

10.00611 


15. 01M1U 

20.0001) 


3.901)0 

_3-»9MU- 


3.90110 

3.903)0- 


25.0000  3.9000 

_ SAMPLE  CALCULATION  !  T«0  INJECTION  b  f  AT  IONS  USEll 


02 

-OjJJL 


II  .05 


-XlulL 


0.000122 

0.0005 


75. 

MiU. 


164. 

0.1)005 


453 .25 


1330.0 


75. 

600. 


0. 00005  0.0001 


J- 


SAMPLE  CALCULATION 


For  purposes  of  illustration,  selected  pages  of  printout  for  a  sample  calculation  are 
presented  below.  The  duct  geometry  is  given  in  Fig.  10,  and  two  injection  stations  are 
considered.  The  initial  conditions  can  be  deduced  from  the  printout,  and  the  initial 
integration  step  size  was  0.0001. 


BEGIN  STM  i  Vi  1  RUN  lNf8ftM»T|8N  SAMPLE  CAlCi)l*T|8N  I  TM8  INjECTJBN 

INITIAL  AREA-*  i 18.501 

_  ,08000b-"  "  HG*  .176002  POi  1359.11  T8«  <54.0*9 


A/A0-  1.00000  W V/WNC  •  .12201 

CVT  '.122000-03  VG-  184,000  "  TG."'  '  453.250 

TTi  "T5oa'coo’-'b3"  vl-  ts.aooo  "  ts;  600.000 


rLP*  -.737856-02  VLP- 

'CVpi  ~  1  7  3767  6-02  VGP-  89.70*1  fSPi  2.9*443  PP>  *29,0501 

Xi  .ioobno-02  tiff*  .176101  P8i  l359,l*  T8-  454.090 


•  1330.00 

•  .500000-94  0ELT--144.75 


A/A! 

HHHHH 

m 

IHiBH 

C  V. 

.122714-0 

m 

G» 

1S4.0O0 

TG» 

453.250 

EL- 

.499286-0 

3  V 

L- 

75 . 36  3* 

TS. 

598.379 

1330,08 

.  499762*04  DELT..145.12 


'.6B9IJ9-02 —  VGP- - *9.4455  - TCP."  "2, 77425" - PP. - -2*. 9449 

.200(100-03''-  MG-  .176109""  '  ~  POiT  "  1359.11  T0-"  456.090 


a 

■  .vvyyii 

kx/WNC  —  .IZJ 

JV/-00 

CT« 

•  l Z  3 J  ^2*03 

VG- 

1E4.01T' 

tG» 

453 .255 

EL- 

.495418-03' 

9L- 

75.7326 

TS.  ' 

59  6  r*50 

rvp- 

- . 6  4  6Vy6- 32 

't646B36-02 

VLP- 

VGP- 

3759775 

*9.2535"  " 

TTP- 

TGP. 

-l-BJV. 

2.5*74 

JT- 

~.  300000-03 

HG- 

“  .17611* 

P8. 

1359,11 ' 

17X5" 

CV. 

•  ",999668 

"  .124010-03 

W 

VG- 

V/WNC  . - 7174 

1*4.027 

025-03 - 

TG- 

453.251 

CL- 

.497990-03 

VL- 

76.1054 

TS- 

595,403 

’»  1329, TT  "  ~ 

3.  "  ,'499539-04— DELT. -143,60 


■  1329.99 

«  ,499329*04  DELT--142.15 


.609140-02  VGP- 
.  400000-63  MG- 


*9.0642  TOP."  2.40915 

,"176124  PB-  1359,1*  " 


PP-  -2*. *252 
Tflr  456,091 
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BEGIN 

~r. 

*7*o* 

ST.T13N  2 

.503000-01 

,977904 

RUN  INFORMATION  SAMPLE  CALCULATION 

"MG*  .1*6252  P0."  1359.1* 

WV/HNC  •  • 152594*03 

i  t  T  NO 

TO. 

injection 

456.151 

ev. 

1*8,333 

TUi 

453.17* 

Pi 

1320.61 

~ Tl« 

"•*89420-03 

VL. 

150.347 

TS« 

506.901 

D. 

.409592-04  DELT.-53.725 

~Tlp,— 

■^Tl  a  0506-03 

VLP. 

594.916 

TSP* 

-488.390 

DP. 

-.349416-05 

EL. 

USSS 

75.63*0 

T5. 

600.000 

Da 

.100000-02  CELT. .1*6.02 

t  LP* 

- .  227451-02 

VLP. 

1329.20 

TSP. 

•5148.88 

DP  • 

-.151954-03 

'  cvp.' 

.257909^62 

VGP*  ' 

87.66*3 

tgp. 

-.491097 

PPa 

•20.5993 

X. 

bL^BH 

mini 

pd» 

-rs59:n 

TB* 

45*7203 

A/TToi- 

7956629 

rfV/NNC  ■ 

.178992-03 

— CVi 

.174960-075 

~V6« 

192.701 

TO. 

453.389 

P. 

1327.10 

fL« 

. <66313-03 

vL* 

~ I7D.7T7-" 

rs. 

*90.101 

□  . 

'TLP.  ~ 

— -.37l977-or‘ 

VLP« 

270.1*9 

TSP. 

-229.942 

DP. 

-.129094-05 

fl* 

.*76711-03  - 

VL« 

120.4*7 

TS. 

536.923 

D. 

.984226-04  DELT • -83 . 834 

(  LP. 

-. 13*6 / 6-03 

— TUP* — 

563 .0»6 

TSP. 

-*6o;tv5 

or* 

-  1  V«:60»/-U3 

~"~;i7m3-oj"" 

VGP. 

87.1*38 

*  TSP. 

-2.1729* 

pp. 

-28.8080 

-n - 

m’oimj  — 

'  • 

.188627 

"  P8» 

1359.14 

TO« 

436.228 
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STATION 

•  2 

PUN 

Information 

SAMPLE  C*LCULAT)ON  1 

TWO  INJECTION 

A/AO" 

'7936343 

MV/WNC  • 

,lB4743~ffy 

cvr" 

.161769-03 

V6" 

1977042 

- TG. 

452.973 

p. 

1325.73 

fL» 

VL" 

160.760 

TS" 

“Til  .050 

D« 

.466016-04  OfcLT.-2B.077 

~flp* 

-  •  220*472-01 

VLP" 

170.755 

TSP" 

-148.711 

DP" 

-.771*45-06 

tlv - 

~747237o-03'  ' 

VL" 

141.641" 

TS"' 

520.547 

c. 

.961230-04  DELT.-67.574 

FLP" 

-.5B4431-04 

- 330.517  ' 

iliigHiifll 

-2*0, 957 

CP. 

-,*04670-05 

“eVpi 

.604506-64 

VG  Pi 

87.69*1 

tGP" 

-2.49192 

PPu 

-29.6220 

X* 

.?0i551T 

NG" 

~  7192963 

PB" 

1359.13 

TO. 

456.246 

cvi 

. 167928-03 

'VG*.  ~ 

201 .562 

TGb 

452.840 

P. 

1324.19 

fL*'  " 

T*63»5$-f3r  — 

VL"  " 

1*6.265 

TS" 

474.635 

D. 

.467669-04  DELT.-21.796 

FLP. 

-.l5l?8l-&4 

t$p« 

•lfl5 .883 

DP. 

-.531801-06 

tl. - 

TXTooSj-'OS  _ 

VL"  '  ' 

1557701 

TS. 

510.570  ' 

D. 

.979644-04  DELT.-57.730 

^LP« 

-  .'345209-04“  ‘ 

'  VLP" 

230. *44 

TSP. 

-159.695 

DPi 

-.239803-05 

7496903-04 

smrnr 

HBaum 

-2.69889 

PP" 

Y. 

.252500 

KG. - 

'  .197T79 

P8." 

1359713' 

TO" 

456.260 

A/AO" 

.896670 

HV/MNC  • 

.196096-03 

....  _ 

CV. 

.  19004(1-03 

Td. 

P. 

1322.60 

tli - 

'  .463295-03 

VL" 

194.251 

TS. 

469.925 

D. 

.*87*53-04  DELT. -17.228 

"FLP. 

-.114927-04 

VLP. 

io* . 0*6 

TSP. 

-79.9794 

DP. 

-.403066-06 

FL» 

.966524-03  ’ 

mam 

166.072 

TS. 

503.533 

D" 

'  FLP. 

-.236416-04 

VLP" 

179.221 

TSP. 

-119.060 

DP. 

-.16*709-05 

C'VpV  ‘ 

~  .351410-04 

VGP* 

89.4507 

tGp. 

-2.86146 

PP. 

-31.4426 

X  ■ 

.50)53! 

po« 

1359.13 

TO" 

•27Q 

A/AO" 

.“879049“  ’ 

NV/MNC  • 

.191560-03 

”cv» 

.1915*4-63' 

VG" 

210.433 

TG. 

452.556 

p. 

1321.07 

Fl" 

.462796^07 

VL. 

199.166 

TS. 

466.504 

D" 
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STATION 

■  J 

bun 

Information 

sample  calculation  i 

TWO  INJECTION 

ri_p» 

-.933774-03 

VLP" 

98.5421" 

TSP. 

-63.3097 

DP.  -.320423-06 

fi.r  " 

.467657-03 

VL* 

173.931 

TS« 

490.412 

Da  .977941-04  DELt.-43.036 

flp. 

-.179574-04 

VLP* 

-93*6063 

56 «  =7124969-33 

~rvp« 

.272740-04 

VGP" 

9H.1414 

TOP. 

-2.99409 

PPi  -32.3136 

x« 

73*1700 

MG" 

.206010 

Pfli" 

1339.12 

TO*  456.270 

A/AO* 

706090$ 

r-  •  .  3»» 

cvi 

"7192^92-03  " 

VQ« 

21  r. 066 

TG.~" 

432,400 

P.  l3£9 , 39 

fl« 

.  462340=03  — 

Vl« 

2047081 

. TSi 

463.394 

D>  .407123=04  '  DELT.-ll .194 

FLP" 

:  tffSijKi 

HU1 

- 93  7  2  S  37 

t5p» 

•50 • 5175 

GPi  -.274892-06 

_FL. 

.  466025-03 

"VL*  ~ 

1017693  " 

TSi 

493.954 

D.  .977373-04  DElt.-41.3S4 

FLP" 

-.14047  3-0= 

~VIP« 

130.137 

TSP. 

-79.2413 

DP.  - , 96i042-06 

rvp. 

iiuittum 

56.6390“ 

tbP. 

=3.12331 

PPi  "3J.Z455 

X.  •  — 

.402900 

MG"  “ 

.716319 

- P6.  . 

1339.12 

—  T94  436.200  - 

~ A / AO" 

.043421  ‘  ~ 

WV/NNC '= 

.193031-63 

..  _ 

'  '  . .  . 

cv. 

7193/94-03 

-9Ui 

219. 734 

432 .217 

Pi  1 J  1  7 . 6  / 

-rc, - 

"7*01959-03  ~ 

VL-  - 

~  203.777  - 

- T5i - 

461.243— 

"  C.  .406994-04  DELT>  -9.006 

TCP* 

-.675770-05  " 

VLP  « 

9073494 

- TSP* 

-4173695 

DP  4 - -.237276-06 

re? - 

.  *661V6-U.i 

"TU* 

1377*33 

TS» 

4?D  •  4?®l 

D»  ,  9  7BVJB-U*  0tLT«*4l , UZ* 

FLP- 

-.106320-04 

VLP" 

117.739 

TSP. 

-43.3013 

DPi  -.742663-06 

cvp« 

.173781-04 

VGP. 

91.4931 

TGP. 

-3.24423 

PP.  -34.1703 

xa 

.454100 

MG* 

SUMS 

- P9i - 

1339.12 

T3*  436,292 

A/A0« 

.026362 

HV/UNC  ■ 

.194634-03 

.  ..  -  . 

cv« 

.194616-03 

VG« 

224,434 

TG. 

452.066 

P.  I3i3,09 

FI. 

.461664-1)3 

Vl« 

- 2137377 

TS" 

43V. 320 

D=  , 486880-04  KELT*  -A. 2^2 

FLP" 

-.396296-03” 

VLP"' 

89.2937 

TSP. 

•34.1336 

DP.  -.209622-06 

FL- 

.465697703 

~Vl.~ 

193,220 

TS.  " 

407.134 

D.  .976507-04  DELT.-39.06i 

FLP» 

VLP* 

100 .81 B 

rSP. 

-36.W62 

DPi  -.624233-06 
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station 

*  2 

RUN 

] NF  ORM A  T I  AN 

x« 

1.01090 

MG* 

.26*155 

PO* 

1359.09 

TO* 

*56.323 

A/*0* 

.677729 

HV/WNC  « 

.199661-03 

cv. 

.199*21-03 

KiS 

T  Ga 

*49.56* 

Pa 

1293.61 

fl» 

,  *5  9422*03 

VL* 

2*3.398 

TS. 

*  5o . 1 96 

D* 

.*86091-0*  DELT*  -.312 

fl  p» 

-.207097-05 

VLP* 

9l.5*25 

TSP* 

-7.615*5 

DP. 

-.10*781-06 

FL  * 

.**2932-03 

VL* 

2**. 8*1 

T  S  > 

467.916 

Da 

,97*651-6*  DELT*-1B.o3* 

TLP* 

-.291739-05 

VLP* 

87.4*30 

TSP. 

-21.5639 

OPa 

-.20*776-06 

cvp* 

.586651-05 

VOPi 

9ft. 1360 

tGp* 

-*.**170 

ppi 

-44,6817 

x« 

1  •  26690 

MG* 

,289881 

1359.08 

TO* 

*56.329 

a/a'o*  - 

.6271*6  ' 

HV/WNC  > 

,200997-03 

cv. 

.200956-03" 

VG* 

301. *07 

TG. 

448.709 

Pa 

1281.65 

fl* 

7*58726-03 

VL* 

286,80* 

TS* 

**8.520 

D« 

.465645-04  BELT*  .189 

FL  P* 

-.250109-05 

VLP* 

9o.«648 

TSP* 

-5.59*31 

DP* 

-.662965-07 

Fl« 

. **2293-03 

VL* 

2*7,186 

TS* 

*63.170 

D> 

.97*202-0*  DELT*-1«.*61 

Flp* 

- .21*1*1-05 

&6.46M) 

TSP* 

-15,9636 

DP« 

-.156*21-06 

cvp. 

.*5*063-05 

VGP* 

89.40*2 

TCP. 

-*.5136* 

pp. 

-45.0061 

1.52290 

MG* 

.312016 

P6* 

1359.07 

TO* 

*56,333 

A/AO* 

W V/WnC  . 

.20ZU/2-IIJ 

CV. 

. 20  7o3l-o3 

VG* 

32*. 001 

TG. 

**7,528 

Pa 

1270.12 

Fl* 

.*55132-03 

VL* 

309.157 

TS* 

**7,23* 

.485635-0*  DELT*  .293 

TUpl 

-.215957-05 

VLP* 

B573S5V 

TSP* 

-*.61066 

DP. 

- f  763071-07 

fl. 

.**1811-03 

VL* 

288,777 

TS* 

*59.560 

C. 

,973864-0*  DELT*-12.053 

FLP* 

- « 1*5*88-05 

VLP* 

8?.fr8>3 

TSP* 

-12.3639 

DP* 

-.116608-06 

CVP* 

VGP* 

88,12/3 

TGP* 

-*.77667 

PPa 

-47 .2696 

Xi 

1.75330 

MG* 

.333370 

PC* 

1359.05 

TO« 

*56,337 

A/AO* 

.55*131 

wv/UNC  . 

.202691-03 

Wi  ,?02850-03  V5!  3*5 , 7?B  TTT*  Pi  1258, 11 
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sTat ion 

•  2 

TSFS Rma  T I S  n 

5AH°LE  CALCULATION  t 

FLP* 

"'-.364096-0  6 

VLP" 

iB';r205 

‘  tSP" 

-2.21737 

DP* 

-.200271-07 

FL" 

— .459699-03 

VL" 

4Bl7B9l 

TS  « 

441, 604 

D. 

.472376-04  BELT"  -6.731 

FLP" 

VIP» 

T  rsR" 

-3,30110 

DP- 

• , 231763-07 

CVP" 

.892428-06 

VQP« 

23.0196 

TOP" 

-1.95110 

PP  i 

-17.9706 

Xa 

“*75l810 

MG* 

.495450 

PO" 

1359.07 

TO* 

496,369 

A/AO" 

.402919 

UV/bNC  . 7 

CV« 

”T2o761B-03 

V  G  ■ 

'506.530 

TG"' 

434,793 

P- 

11*8,65 

FL" 

7454735-03 

"VL" 

"  "5037734 

TS" 

435.166 

D" 

.484*32-0*  DELT"  -.375 

FLP* 

VLP" 

9.14606 

TSP- 

-i.5$SS5 

DP* 

-.178683-07 

FL" 

.459619-03 

Ti“"  ~ 

488.293 

TS" 

440,617 

D> 

.972320-04  DELT*  -6.024 

~"fLP* 

-.289989-06 

I. 

•» 

1 

20.3574 

tsp* 

-2.60500 

DP" 

-.20**90-07 

■HUSH 

.792849-06 

VGP. 

-8,86382 

TEA. 

.387965 

PP- 

5.35070 

X  B 

'4.7741(5 

MG*  ~ 

.'494926 

50. 

1359.07 

TO" 

*56.370 

A/AO"  “ 

~  *  7403233  ' 

WV/MNC'T  , 

20  783^2-03 

cv» 

. Z07809-03 

VO" 

508.424 

TG- 

434 , 637 

P  ■ 

114V  ,04 

FL" 

~.4yi{ij-07 

“VC" — 

- 50r.l2T~ 

1 

-J 

w» 

■ 

434.632 

D. 

.*8*369-04  DELT"  .004 

rtp* 

-,4V57jr-7l« — 

~  VLP«“ 

—  3751890 

TSP. 

-1.10352 

DP. 

-.161159-07 

— res - 

.  *5V5«9-0 $ 

VL* 

- 4927555 

1 

4*0.159 

0- 

.9/22/1-04  DfcLT"  -5.3ZJ 

'  rLi**1 

-7237961-06 

■  VLP« 

13.9148 

TSP" 

-2.37062 

DP. 

-.181923-07 

“  CVP"~ 

.713392-06' 

VGP« 

.862684. 

■02  TOP. 

,199743 

-02  PP. 

-.129745-01 

X* 

3.00450 

. 494925 

Pfl" 

1357.06 

TO* 

456,370 

—  A/AO" 

~  .403233 

UV/UNC  ■ 

,208011-93 

evi 

'  ‘.207968-03 

VG" 

506.426 

TG. 

*34,637 

P- 

1149,04 

FL" 

.454511-D3 

VL" 

505. 724 

T5- 

134 .609 

D- 

, 184353-04  DELT ■  .228 

“  FiP* 

-.432021-06 

vup* 

1.85835 

TSP" 

-.644245 

DP* 

-.153462-07 

“TL. 

.459492-05“” 

Vl» 

495.345 

TS. 

439,647 

D* 

.972231-04  DELT"  -4,810 

Flp« 

-.236772-06 

10.4983 

-2,08465 

DP- 

-.166994-07 
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